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The FORTRAN IV computer listings for each of the radar 
cross section models are located at the end of each of the 
appropriate sections L.1 through L.7. Although the programs 
have not been integrated into the Radar Simulation Model, 
they can be easily made to interface with RADSIM. In the 
interim they can be run as independent FORTRAN subroutines, 
which are stored on disc file on the RADC HIS-635 computer 
under USERID-CLEARY and account number 017073100380 in the 
following files:” 


SECTION NAME SUB-FILE LINES 
L.1 Cylinder with Spheroidal RCSM3 10760, 11430 
Caps 

L.2 Hemisphere-Cylinder RCSM3 11440, 13680 
L.3 Cone, Cylinder or Frustrum RCSM3 13690, 17460 
L.4 Thin Wire RCSM1 50, 3120 

L.5 Frustrum-Cylinder-Frus trum RCSM1_ . 3130, 7580 
L.6 Cylinder-Frustrum Combination RCSM1 7590, 12170 
L.7 Stepped Cylinder RCSM2 50, 3710 


The outputs for all programs consist of four linear 
arrays - EVVR, EVVI, EHHR and EHHI - which contain the real 
and imaginary parts of the vertically and horizontally pol- 
arized back scattered fields as a function of frequency. 

The fields are computed at discrete frequency increments 

(DF) in megahertz, between user defined limits. The low and 
high frequency limits are equal to (NMIN-1)*DF & (NMAX-1)*DF res- 
pectively, where NMIN and NMAX are dimensionless indices. 

The magnitudes of the arrays - EVVR, EVVI, EHHR and EHHI - 

have been scaled in order that the target cross section (in 
Square meters) can be obtained from the sum of the squares 

of the corresponding real and imaginary values of the array. 


A typical use of SUBROUTINE TARGET is shown in Volume i 
Part 3, Section 8.24, All subroutines in this appendix can be 
run using this same FORTRAN program, labelled RCSSP. 


L-1 


L.1 CYLINDER WITH SPHEROIDAL CAPS (RADCAT) 


The far-field scattering from the RADCAT vehicle, 
a large cylinder with spheroidal end caps, has been formula- 
ted using an expression involving a physical optics approxi- 
mation of the backscattering from a cylindrical section and 
a geometrical optics approximation to the scattering from a 
spheroid (Ref 1). The resulting expression of the scatter- 
ing field is the following: 


e(@) = 2 \ kasin@ {Seats icose he i(2kasin@ + 7/4) 


gph , L ae 
a | . Tac oan sin Ota" cus 6) 
a sin 80+ c cos 9 


where the physical dimensions of the target and target 
geometry are shown in Figure L.1-1 and k = wave number = 2z7/). 


The computer program subroutine is used to compute the 
field backscattered from the target for the case of verti- 
cal and horizontal polarizations with respect to the target 
rotational plane (defined by the target axis of symmetry and 
the radar line of sight). It should be noted that for the 
case of this particular target, the formulation of the 
scattered field is polarized insensitive. The inputs, out- 
puts, restrictions, and definition of key terms of this sub- 
routine are presented in the following paragraphs. 


L.1.1 Inputs 
The subroutine inputs are read from cards or passed in 
a common block, The used inputs passed in the common block 


are the following: 


~~ © Carrier Frequency in GHz, not necessarily the 


mid frequency 


DF = Frequency increment in MHz 


L-la 
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Fig. L.1l-l RADCAT TARGET GEOMETRY 


L-2 
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NMIN Minimum frequency index, i.e. 


NMAX 


Maximum frequency index, i.e. f,,, = DF*(NMAX-1) 


The inputs which are read from a card are the following: 


MATH 
SYMBOL VARIABLE DEFINITION COLUMNS 
Card 1 a XA Radius of cylindri- 1-6 
cal section (feet) 
c XC Depth of Spheroidal 7-12 
Cap (feet) 
f XL Length of cylindri- 13-18 
section (feet) 
ts) THETA Azimuth angle 19-24 


(degrees) 


L.1.2 Outputs 


The outputs consist of four linear arrays EVVR, EVVI, 
EHHR, EHHI which contain the real and imaginary parts of 
the vertically and horizontally polarized back scattered 
fields (in meters) at frequency increments of DF between 
DF*(NMIN-1) and DF*(NMAX-1). The bandwidth of the output 
then equals DF*(NMAX-NMIN). 


L.1.3 Restrictions 


L.1.3.1 Physical Dimensions 


As high-frequency approximations were utilized, all 
physical target dimensions should be large with respect to 
the wavelength of the illuminating field. Thus, XA, XC, 
and XL should be larger than the wavelength of the smallest 
frequency at which a computation is obtained. The condition 
XC< XA is also implied in the formulation of the problem.” 


L.1.3.2 Output 


The output arrays are passed in the argument list and 
a value is computed only for array locations from NMIN to 
NMAX . 


*also, 1S NMIN < NMAX < 512. 
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DF*(NMIN-1) 


L.1.3.3 Azimuth 


The azimuth angle is restricted to the region between 
0 and 90 degrees. 
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TERM 2 


TERM 4 
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L.1.4 Definition of Key Terms 
Used in Subroutine 
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L.2 HEMISPHERE - CYLINDER 


The far-field scattering from a hemisphere - cylinder 
has been formulated using the Ruck-Ufimtsev formulation of 
the scattering from a cylinder and a modified expression 
has been utilized in describing the hemispherical returns 
(Ref 2). The resulting expression of the scattered field 
is the following: 


iy = a_.i2kh cos@ | 2 . 
had ideale | Sirs E Big. © t___\Bo1- 
0.5+cos— 
3 


5 1 4-3] 
+ Q ‘2 Bs,. + (-———— IS: BB 
a oli 0. 5+cosat7=28) o1+} 


3 to 
z : tan 8 Bo,~e i2kh cos@ 


CW 
*Os(Sy 4. *) 
O 
where bole - JQ) ae by tg 
Byte. Med ai ah 1/9) 
& = 2 Ln wd 
Qe = 20%ks(O - />) 
Q, = QC 7a crs, =< G5) 
k = 2-/X = wave number. 


In the primary equation, the upper and lower signs designate 
the vertically and horizontally polarized returns respectively. 
The factors S and CW are the sphere specular and spherical 
creeping-wave returns respectively. They are expanded as 
follows: 


L-9 


a i -12k(ath cos@) 
S = -_— 
2 (2 2ka : 


a 0.5455 194489 | -4 ; 
Cw, = $(2.7162)K | (1+ 225929) 4 4 O-PU5BP | .-i2kch cose 


exp { -(2.20075K-2-42225 ) +4 (rkatZt1.27065K+223207_ | 


where K = (kay! /3 and the geometry of the problem is shown 
in Figure L.2-1. The inputs, outputs, restrictions, and def- 
inition ot key terms in the subroutine are presented in the 
following paragraphs. 


L.2.1 INPUTS 
The subroutine inputs are read from cards or passed in 
a common block. The parameters passed in the common block 


are FC, DF, NMIN, NMAX and are described in Appendix L.1.1. 


The inputs read from a card are the following: 


MATHEMATICAL 
SYMBOL VARIABLE DEFINITION COLUMNS 
Card 1 a A Radius of hemis- 1-10 
phere (inches) 
h H 1/2 length of 11-20 
cylinder (inches) 
] ASPECT Azimuth angle 21-30 


(degrees) 


L.2.2 OUTPUTS 


The data base output consists of four linear arrays, 
EVVR, EVVI, EHHR, EHHI, which contain the real and imaginary 
parts of the vertically and horizontally polarized back- 
scattered fields (in meters) at frequency increments between 
NMIN and NMAX, 


Fig. L.2.1 HEMISPHERE=CYLINDER TARGET GEOMETRY 


In addition to the data-—base output, if the print option 
is selected, the frequency, and cross-section (in dBsm) and 
phase for the case of both vertical and horizontal polari- 
zation will be listed. 


L.2.3 Restrictions 


L.2.3.1 Physical Dimensions 


All target dimensions should be large with respect to 
the wavelength of the illuminating field. Thus A and H 
should be larger than the wavelength of the smallest fre- 
quency at which a computation is obtained. 


L.2.3.2 Output 


The output arrays are passed in the argument list and a 
value is computed and stored only in locations NMIN to NMAX. 


L.2.3.3 Azimuth 


The azimuth angle is restricted to the region from 0 to 
180 degrees. In addition, an azimuth of 90 or 180 degrees 
should not be used as there are numerical problems in com- 
puting the specular return. If the magnitude of the specular 
returns are desired, azimuths of 89.99, 90.01 or 179.99 de- 
grees should be used, 


L.2.4 Definition of Key Terms Used in Subroutine 
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ARGMNT*PHASE 3 = {1 ae Aelia cos@) 
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0.44525 
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PART 1 = 8(,7162)K exp { -(2.20075K- ) 


PHASE 4 = exp {+4 ( mkatFtl.27065K+2:25207 } « g -i2kh cos@ 


0.54555 


0.94489 
YY = ii ey ee 


Q PLUS = Q4 = Q(2ka(@ - 7/2) 


QMINUS = Q_ = Q(2ka(7/2 - @6)) 
L.2.5 Subroutines Utilized 
Subfunctions: 


Q(X) computes the exponential smoothing function Q 
for real argument x 


Subroutines: 


BESL (TERM 2, 
JO (TERM 2) 
J1 (TERM 2) 
J1 (TERM 2) 


J2 (TERM 2) 


XJO, XJ1, XJ2) returns 
in XJO 
in XJO 
in XJl 


in XJ2 


| 


SUBROUTINE TARGE! CEVVR, EVI, EHHR, EHHI9) 
+ * TARGET ST-S, HEMISPHERE CYLINDER, UFIMTSEY TECHNIGUE * »* 
COMMON MOVER. M. NMIN,. NMAN, OF, FC. Fa, T@ 
NMI? MINIMUM FREGUENCY SAMPLE 
NMA = MANIMUM FREQUENCY SAMPLE 
CF FRECUENCY INCREMENT IN MHS 
FC = CARRIER FREGUENCY IN GHZ 
DIMENSION EVVRCS425, EVVICS412)9, EMHHRS S429, ERHI (S12) 
FREG@M S412), SIGMAVCS12), SIGMAHCS12), FPHASEY(812), FPHASEHCS12) 
MELEX COEFFI, COEFFe, FHASEL, FHASES,. FRASES, PHASE4, FFYYV, FFHH 
FEWW4, FFHHL FRYVe, FFHHe, FRFVYE, FFHH:, FFY¥4, CFYY, CFHH 
q = RADIUS CF HEMISPHERE 
ri = HALF THE CYLINGER LENGTH 
ASPECT = AZIMUTH ANGLE 
itd FRINT GETICQN 
t+ DIMENSIONS ARE IN INCHES AND ANGLE 15 IN DEGREES #* * 
READS, 16603 KD ASPECT, M4 
FORMAT Sd Se ae rae 
WRITE ¢€6. 44414) Ay H, ASPECT, M41 
FORMATCAHG, Fi ety Pas Se Say Hee 4 Fos. S é¢ 
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C = 11. &676Sa76 
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COMPUTE @ CSMOCTHING) FUNCTIONS 
ZF = KEKCA * (THETA - PICGVRe) 
ZN = K2KCAH * CPIOYRe - THETA) 
OPLUS = OCZP) wwe 
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FFVYe = COEFF2 # PHASES 
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HEMISPHERE SPECULA® 
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HEMISPHERE CREEPING WAVE 
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FUNCTION C(Z) 


IF ( 2-GT.e- 2-) GO TO 10¢ 
IF € ZeLTe—2e} GO TO 20 


hé = ABS(Z) 


9 
" 


LeG/(1.60 # 647047#AZ) 

Y = 160 = P#(.3460262 — PH .LOG587496 — .747TESSE*P))*EXP(-AZ AZ) 
IF (2) Ze46 

© = (166 = YD/26 

xe TURN 

& 2 65 

SE TUKN 

& = (1.0  VY)/2. 

RETURN 


© = le 


= Ue 
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04/26/76 AUTUFLOW CHART SET - FWO/SCL RADSIM 


CHART FITLE — NON-PROCELUKAL STATEMENTS 


COMMON MOVER, My NMINy NMAXy DFy FC, PW, TC 

. DIMENSION EVVR(S12)5 EVVI(S12), FHHR(S12), EHHI(512) 
» FREQ(5S12)— SIGMAV(®12), SICMAH( 512), PHASEV(512), PHASEH(512) 
COMPLEX CCEFFL, COEFF2, PHASEL, PHASE2, PHASE3, PHASES, FFVV, FFRH 


9 FFVV1», FFHH1, FFVV2, FFHH2, FFVV3_ (FFHH3, FFVV4, CFVVy CFHH 


’ YY 
1000 FORMAT( 3F10.0,15) 
1111 FORMAT(IHC, © A = %, FI565 5 ICXe 8 H = Sy FISS re 


1HOy * THETA = %) FlSe5y 15X_ " PRINT OPTION = *, 15) 


L.3 CONE, CYLINDER OR FRUSTRUM 


A generalized expression of the far-field scattering 
from a cone, cylinder, or frustrum has been formulated using 
the Ruck-Ufimtsev technique (Ref. 3). The resulting ex- 
pression of the scattered field is the following: 


{-ay [ [Jo (sq )+45, (>) C(n_)+[5,(6,)-45, (6 ,)] BOn_, +9) ] QC) 
-ay [ [426 6,)-15,(£,)] (n+ [59 ( <4) +43) (61)] B(n_,2-8)] Q(4) 
cay [ [22¢é2)+43,(49)] : oe 
2 | [32 <o)+45, (52)] C(my)+ [59 (Zp) 15, (Ey) ] B(ny,+0)] 
~ay [ [5p (54-454 (6 )] c(n,)+ (5,65) +45, (Z5)] B(n,,2-0)] 2” 
2 LLJ20$2)-45, (¢5 HIE ol bg) 425) (59)) Bin, ,a-8)]e° Q(3) 
+a J X€. pat J. (0) dean aw+Q : 
ta, [ 5, (6,) 445, (5, )] Fean(a-0) FC +, )Q(34) 
ta, [5,(4,)-44, (¢,)] dean(otoyr(+,yeMQc1) 


+a, [ 5,4,)+49, 6G, )] Leana-0)F(>,)eiQ¢34) } 


and for 7/2<6@<7 , 


{-ay [s2¢ 61) +43, (<)] C(n_)+|5,(¢4)-15,(¢,)] B(n_,7-a-0)]Q(1) 
-ay[ [39 G9) +iJ, (5p) ]C(n,)+ [5 ( fy) -i5y (£y)| a(n, ,32-@) Jet 
8g [ [¥p(%p)-£5, (4) (ny) +[J, (69) 443, (£5)] 660 9-2) Je’ "Qc3) 
+a, | Joby) -45, (41) | Stan (a+9)F (71 Q(1) 

Fa, | Jol 42) -id, (4p) | dean (+o) F(r9)e"Q(1) | 


where the geometry of the problem is shown in Figure L.3-1 and 


C(n) = : sin “eos = -1 ee 
B(n,¢) = £ sin Sleos = -cos 201 
9 n nt n n > 
74, = 2ka,,,csca cos (a+ 9), 
(>) ee 
ae 2ka cscacos (a@— Q@) 
(7) 2 ’ 
4 1 
¢, = 2ka; sin 0, 
% = 2kh cos 8, 
eee 
"+ 2 = ea 


Q(1) = Q(2kay (7-0-0), 


k = 2n7/X = wave number, 


—— Ne a NALA eS Aca ee 
m lS 


o = 
o 
n 


Fig. L.3-1 CONE, CYLINDER, OR FRUSTRUM SCATTERING GEOMETRY 
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Pa he 


Q(3) = Q(2kay(a-8)(9-77/2)) , 
Q(34) = Q(2ka,(«-0)), 


= Bs 
Q¢4) Q(2ka, (5 @)). 


The above solution to the scattering from a frustrum can 
simply be reduced to that of a right circular cylinder by 
setting F(Ti) = 0 (i.e., @ = 0). For the case of a large 
right circular cone, it is important to note that numerical 
difficulty will be encountered if the above expressions were 
used to compute the scattered field near and at the conical 
specular aspect, @ = 7/2 -@. An asymptotic solution to 
this problem can be expressed as 


; = ip 

lim {Voe'#} =+Vr ase ; 
@— 7/2-a Vv 
H 


{- [ 39(¢y)#45 (Lp) JoCay dt [F,(4y)-45, (65) | * 


[iSkazescasin(a+e)}} : 


The inputs, outputs, restrictions, and definition of 
key terms in the subroutine are presented in the following 
paragraphs. 


L.3.1 Inputs 


The subroutine inputs are read from cards or passed in 
a common block. The parameters passed in the common block 
include NMIN = minimum frequency number, NMAX = maximum 
frequency number, DF = frequency increment (in MHz) and 
FC = carrier frequency (in GHz). The inputs read from a 
card are the following: 


L-27 


MATHEMATICAL 
SYMBOL VARIABLE 


Card 1 ay Al 
a9 A2 
ho H2 

KONFIG 
(e) THETAD 


print option KP is positive. 


L.3.3.1 Physical Dimensions 


DEFINITION 


Smaller radius of 
frustrum (inches) 


Larger radius of 
frustrum (inches) 


Height of frustrum 
(inches) 


=]1 Frustrum 
=2 Cone 
=3 Cylinder 


Azimuth angle 
(degrees) 


L.3.2 Outputs 


L.3.3 Restrictions 


COLUMNS 


1-10 


11-20 


21-30 


31-35 


36-45 


The data base output consists of four linear arrays, 
EVVR, EVVI, EHHR, EHHI, which contain the real and imaginary 
parts of the vertically and horizontally polarized back 
scattered fields (in meters) at frequency increments 
spaced DF MHz from NMIN*DF to NMAX*DF, 


The frequency, vertically polarized cross section, and 
horizontally polarized cross section are printed out if the 


All target dimensions should be large with respect to 
the wavelength of the illuminating field. 


L.3.3.2 Output 


The output arrays are passed in the Argument List and a 
value is computed and stored only in locations NMIN to NMAX, 


L.3.3.3 Azimuth 


The azimuth angle is restricted to the region from 0 to 
180 degrees. In addition, although the formulation is mathe- 
matically valid at specular points, the correct computational 
results are not provided at these angles. Thus a small angular 
offset should be used at angles of 0, 180, and ( 7/2-a) degrees. 


L.3.4 Definition of Sample Terms 
Used in Subroutine 


FSTCOX = C(n) = 2 sin = {cos a ae 
n n n 
we ee 
Where m,, 7] z = 
FSTCON = C(n_) 
FSTCOP (n,) 


=C 
TAUSQ |}t- mad ¢ CSC 4 
J 


XOPX11 = 


XOPKX12 = J,(f9) HJ} (5) 
Qi = QC1) = Q(2ka, (7-2-0 ) 


SECND 1 = B(n_, 5 + @) 
WHERE = B(n,@) = 3 sin 2 
n n 


-cos 26 


n -1 
Tr n 


cos 


V 
TERM 1 Ii 


if 
--a) [ | 5265p) +45 (4) CCD [9,0 64)-45 (£)] BCn_, F409] (2) 


TERM 4 tea) [ [42661-8592 (61) 669+ [Fo Cp) +451 1)] BC, 2-0)] Q(4) 


rer 2 |¥t=-a2 | [32 éod+433 (623) o¢04% [96 ¢2)~19 (49)] 8g 040)] 
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TERM 3 Ee nay [ [Jp(bp)- 15, Ep) ]o(ny r+ [Fy (p) +45 (29)] 
B(n, ,a-@) eo" Q(3) 
TERM 5 {= +a, [ 1,(¢,)- 49, (Gd) ean rorC a 


WHERE = the upper sign is used for V polarization and 
the lower sign is used for H polarization 


2 [5 (fp) -45, £9) | {5 [s§kazescasin(a+9) | 


YY =~ [ 34 (Zp) +43, (62)] O(n) 


L.3.5 Subroutines Utilized 


Subfunctions: 
1. FIRS ¢ XN) computes c(n) 
2. SECO (PHI, XN) computes B (n, @) 


3. Q (X) computes the Q function value with 
argument x 


4. F (TAU) computes the F (T) function value 
Subroutines: 


BESL (ARG1, XJO, XJ1, XJ2) returns 


Jo (ARG1) in XJO 
Jl (ARG1) in XJ1 
J2 (ARG1) in XJ2 


- 


oO 


SUBRCGUTINE TARGET « EVYR, EWVVI. EHHR, EHHI, THETAD) 


* £ * &# & GENERALIZED PROGRAM FOR A FRUSTRA, CONE, CF CYLINDER + * # 
IFIMTSEY SOLUTION FOF ti + 


? 
ee = FRONT-ENT + 
A. EF ~END GE At} NHCHES * 
z TOTAL EN 4 
FHE ' F (CEGREES> + 
} PIER FI HENS CGH t 
[ r E 1ty INCREMENT “DEGREES. - €. GE 6.4 + 
KONF TC AFGET CO} SURAT ION + 
rR + 
‘4 + 
= ue INCDER 4 
+ 
* £ 4 ALL CIMENSIONS ARE IM INCHES ANG ANGLES ARE IN GCEGREES to Ft 
+ 
FE Pid: Te 
EHHRF OC S129, EHH este 
SIGMAHL S123 
Pale GQMUXTe, KEPRIL, KEAN Li 
TERMLH, TERMEY, TEENGH, TERNEY, 
TERMSY TERMSH, TERMEY, TERMEGHy 
TERMGSH, FFVY. FFHH, F, FTAUS 
wey EY ng! 
+ 
+ 
t 
REAOCS, 1660) Ai, Ade, He, KONFIG, KI 
LOCH FORMAT SEES 2s. > 
AL = fie Ad 
ALPHA = ATANSCALP. Hed 
FI = Z. 141592652589 r9 
ANGLE = ALPHA * 168.8) ¢ FI 
* 
EF CRONFIG = 23-36, 46, ‘86 
* 


2G WRITE (6, 3030) Hz, fil, AZ, ANGLE 

AZ@ FORMAT (14O0/7° FRUSTRA CUFIMTSEY SOLUTION)’ 4-7 FRUSTR 
1A LENGTH = °,F4@.6,° INCHES’ -’ FRONT-END RADIUS = °,F16.6,° INC 
ZHES’ /’ BACK-END RADIUS = °,F4@.6,° INCHES’ ¢’ FRUSTRA HALF-ANG 

SLE = ’,F10@ 6,° DEGREES’) 

GO TO 6@ 


40 WRITE (6, 2@4@) He, Ae, ANGLE 
364G FORMAT “1H@/ 7’ CONE CUFIMTSEY SOLUTION)’ “-° CONE 
1 LENGTH = ’%,F1@.6,° INCHES’ 7’ BASE RADIUS = °,F1@.6,° INCH 
ZES’ ¢” CONE HALF-ANGLE = °,F1@.6,° DEGREES’ > 
GO TO 6@ 


* 
S@ WRITE ¢6, 3656) He, Ae L-3) 


ow 2656 FORMAT ¢1HG//-'’ CYLINDER CUFIMTSEY SOLUTIOGND* 47% CYLIND 
1ER LENGTH = °,F41@.6,° INCHES’ /* CYLINDER RADIUS = °,F16.6,% IN 
2CHES’ > 
C * 
6@ WRITE ¢6, 795 THETAD 
79 FORMATCLHG, * THETA = * . E15. & tLe 2 
c * 
c = 141. 8@eesa7s 
OTR = PI “¢ 186.6 
PISGRET = SQRTCPRID 
PICVRe = PI ¢ 2. @ 
Al2 = 2@ Gk AL 
Ace = ¢.@* A2 
Hee = 2. @* He 
KNPOS = ¢€2.@ ¢-2. 8) + CALPHA 7 PI) 
KNNEG = (2°. Sos = CACPHA- + PTD 
C COMPUTE CCN? 
FSTCOP = FIRSCHNPOS) 
FSTCON = FIRS¢XNNEG) 
THETA = THETAD * DTR 
STHT = SINC THETA? 
CTHT = COS¢ THETA? 
Le AFT = ALFHA + THETA 
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L.4 THIN WIRE 


The far-field scattering from a thin wire has been pro- 
grammed using the formulation of Ufimtsev (Ref. 4 and 5). 
The expression includes the higher order scattering terms 
which arise from waves which are launched from one end of 
the wire, traverse the wire length, and are diffracted into 
space and reflected along the wire upon reaching the oppo- 
site end. These scattering mechanisms and the geometry of 
the problem are shown in Figure L.4-1. The ray components 
of the scattered field are described in References 4 and 5. 
The expression of the scattered field, which includes only 
the horizontal polarization response, is the following: 


560) 2/7 m 0 _ 
k sin 6 sin 26 [in (<~,)| 
yka sin 6 


where Y = 1.781 


k = 2n/X = wave number 
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Fig. L.4-la GEOMETRY OF A THIN WIRE 
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L.4.1 Inputs 
The subroutine inputs are read from cards or passed in 
< a common block. The parameters passed in a common block 


include: 


NMIN = minimum frequency sample number 


NMAX = maximum frequency sample number 
DF = frequency increment (in MHz) 
FC = center frequency (in GHz), 


The card inputs are the following: 


MATH 
SYMBOL VARIABLE DEFINITION COLUMNS 
Card 1 Q THETA Azimuth angle 1-10 
(degrees) 
L XL Wire length (inches) 11-20 
2a DIA Wire diameter (inches) 21-30 
L-52 
ET ET 
verre 


L.4.2 Outputs 


The data base output consists of two linear arrays, 
ETTR, ETTI, which contain the real and imaginary parts of 
the horizontally polarized backscattered fields (in meters) 
at frequency increments spaced DF MHz from NMIN*DF to 
NMAX*DF., 


If selected by setting KP to 1, the wire radar cross 
section (in dBsm) versus frequency will be printed out. In 
addition the real and imaginary parts of the exponential 
integral, computed in subroutine EXPI, will be printed out 
{fF EP <0, 


L.4.3 Restrictions 


L.4.3.1 Physical Dimensions 


The wire length should be large with respect to the 
largest wavelength of the illuminating field. However, the 
wire radius should be much smaller than the smallest wave- 
length of the illuminating field. 


L.4.3.2 Output 

The output arrays are passed in the argument list and a 
value is computed and stored only in location NMIN to NMAX. 
L.4.3.3 Azimuth 


The azimuth angle is restricted to the region from 0.5 
to 89.5 degrees. 


UW.4.4 Definition of Selected 
Terms Used in Subroutine 
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L.4.5 Subroutines 


The subroutine EXPI (ARGZ, EIXR, EIXI, KP) is used to 
compute the real (EIXR) and imaginary (EIXI) parts of the 
exponential integral? 

rer =f wat 
nd 


where the argument X is passed as ARGZ. 
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6 i6. Sf 
¢ IDENT BECAGDG1. HANCOCK, G17@°3106786,. RCSM1L 
£ GQFTIGN FORTRAN 
$ FORTY LSTIN, XREF, MAP. DECK 
£ LIMITS G@&, 59K, @, Sk 
SUBROUTINE TARGET <CETTR, ETTIL HR. KIT) 
C 
C ** THIN WIRE CW RESFONSE * @.5 TO 69.5 CEGREES ASPECT ANGLE** 
C + SOLUTION BY UFAMTor oO, SIMPLIFIED EY HONG FOR BACKSCATTER + 
C * REFERENCE - SHORT PULSE SCATTERING EY A LONG WIRE, BY & HONG 
C * TEEE ON AP, VOL AP-16,N0. 5. MAY 1968, PP. 28-342 
C 
C 
COMMON MOYER, M.NMIN, NMA. DF. FC, Pa, Th 
C NMIN = MINIMUM FREQUENCY SAMPLE 
C NMAS = MAXIMUM FREQUENCY SAMPLE 
C CF = FREG@UENCY INCREMENT IN MH2 
c FC = CARRIER FREQUENCY IN GHe 
C 
COMPLEX S5¢1aa, 
A ACLO, ALN. CS. ACLZC, ST1, ERA, EXOP. EXOM, EXAP. FSIA 
& EXGPF, EXGMP, FPSIF,PSIM, ACL, FCI, ACe.SCO1, GCL. GC4 
C SCb2, BCL, SCOs. KCI. EC1,.SCb4, SCD, FRE, FRO. FPLG, 
C SPL, SPS. SF) 3F SP, SFL, SSS, RSG VERSU 
E o OR Ps Pidde 
F »CCERL. 16 »SB4 
i) »EXGA, Exo 
1 RCS, OCS 
REAL FROCS125, ETTRCS42>. ETTICSLe). SIGMACSI2>, KR O512), X1CS12) 
REACCS, 16605 THETA, L. CIA, KP 
KF = i, PRINT GUT ROS VERSUS FREQUENCY 
1668 FORMAT ¢3F16 3. Tes 
WRITECE,11G@> THETA, XL, OIA 
4208 FORMAT <6. ASPECT ANGLE = “) FY, 2. WIRE LENGTH 2 .%. FS. . & es 
H WIRE PIAMETER = PRS we Fe ( LENGTHS ARE IN INCHES 
B) 
SF CTHE Ds. Sti sake GR. THETA GT, 82:35 ) Go TO) Seu 
RAC = DAVE & 
FI = 3. 1415926 
PIJ = CMPLKXC@ G, PID 
PIage= Fl/’e. @ 
PIJGe = CMPLe.G &, FIGe) 
GAM = 1 fe1GFre 
TH © THETA ®© (1286.0 ¢ PID 
CT = COS*TH) 
CTGQe = COS¢<TH/e. &) 
CTGeS= CTGe * fTOe 
CTaerF= CToes * CTOes 
ST = SIN «TH? 
STT = SINC THe. @) 
STQe = SINCTH/(e. &> 
ST0eS = .SsT0e * STOe 
STOZF = STO2S * STOZS 
GeCT = 1.64 CT 
OnNCT = £.6= CT 
AKTH = 2 G@/CST#STT) 
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OO €@@ IFW = NMIN, NMA 
KI IFW - fi 


FREQ = XI * OF / 1000.0 L-S$ 
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aAmoOno 


ANAC 


FRC IFW) = FREQ 

BKO= Stei44id + FRED D 

XKLE = &K* KL 

XKA = KK* RAC 

GkKA = GAMN* KA 

Ci = @ @“¢ GKA * ST) 

ACL = ALQG ¢ C4 2 

ACLC = CMPLS ¢€ ACL» PIGS > 

CS =¢¢@ @ 1. @)¢CACLO*ACLO))*¢-AK THD 


++ DETERMINATION OF INTEGRAL AND CTHER TERMS 


ORDER SCATTERING 

CSL = ALOG* GKAD 

Cer = ce B/HKAD*CKKLAGKAD 

C6RPL= ALOGe CéR) 

CCERL = CMPLACCERL. FIQe) 

EAA = KKLe2. G 

CALL EXPY EAA, EXRA, EXIA, @ 9D 
EX@A = CMFLAX* EXRA, EXTAD 

Eic= COS! ERA» 

E1S= SIN« EAA) 

RCL = CMPLR CELC, &15) 

eC = EMPL ELC, <805 9 
FSIA = APIS = 2. G#CSL) COCERL - EXGARC 
EQFPA = MEL *OMCT 

ALL EXPI CEGRA, EXR, ExI. @ » 
EXOP= CMPLHC EXR, EXI 3) 
Eel =COS¢ EOF) 
ees =SIN¢ EGFA)D 
EXOPP = CMPLAC E2l. -E2S) 

CTP = GKA+GRAeOMNCT Ze @ 
C7PL = ALOG CTP) 
POSTF = CPIS ~ CHPLICCCERL - EXGPeEXORP) 
EQMA = XKL#*ORCT 

CALL EXP! EQGMA, EXRM, EXIM, @ 35 
EXOM = CMPLHC EXRM, EXIM) 

fe = COSCEGMAY 
O3 = SINCEGMAD) 

@CL = CMPLHC Ge, OF) 
ocg = CMPLEHC Ge, -OS 3 

C7M = GKAt+GKA*OPC Te © 
C7ML = ALOG<CC?M> 

PSIM = (FIT = C7MLI ¢ CCCERL - EXOMeGCo) 


* UTILIZATION GF FACTORS IM SCATTERED FIELO 


** FIRST ORDER SCATTERING ++ 


C2 = 1.6 GEKA * STOeZS) 
ACLé@ = ALOG (Ce) 

ACL2C = CMPLR ¢ ACLS, FIOGe » 
ST2 = =STO2F * ACLEC 


+ RETURN FROM TRAILING ECGE CF WIRE * 


AL = 1 6/CGKA* CT02S) 

Ae = ALOG*AL) 

ACL = CMPLAXC A2. PIOe) 

Pi = <. Ge KKL*eCT 

Pe = COS¢P1) 

P2 = SINCP1) 

PCie CMPLKC Pe. FF) 

ACe = CMPLHKC CTOZF,@. &> 

SCO1 = ACI*FCI*AC2 
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a = -——5 Io 
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ae 


RMR IO 


USED 


EXPRESSION * 


bp 


IN MULTIPLE ** 


eetindibeendiinemtee deta 


a 


ey 


Pal 37 be) 


eel 


— 
" 


E 1. @/¢GRA#CTOe> 

Re = 2. G@tALOG¢b1S 

BCL = CMPLXC B62, PID 

SChG =C(-COCL#ACSH# Ci eR SIF) 2 


DL = 41. @“C(GKAeSTOS) 

Ce = 2 GeALOUGe 1) 

OCL= CMPLECES, FID 

SE4 = CMPLE*C STOSF, &. &3 
ChE = SBdt #8 GCL OCL HP SIM: 


= 1. @GkKA 
t = ALG E> 
L=. CMP bak Phas > 
Ch4 = CTH#CRCI#ECLePSIPeP SIP) 
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ao 
o 
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‘ 
oo 
ran 
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n 
‘ 
ro 
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f 
O° 
+ 
nm 
mr 
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BS 


FPL = COSC4 GeKKLD 

FF TMCS. Bh 

Fr? COStes GHeHkKL eG CTS 
Cea CINke & RL eee Ty 
PRPS: = EL + e+ | 

FRE = Cie FBS 

FRY = MPR ES 

FRE = CMPL oat ca SP a) 
FPS © Chir PRs. FRA 

FEI@ = CMF FP Ga \F Py 


Fe POD PAP SD Per PS 


‘ 
” 
rh 


SF = (SFL ESF oer Se CTH Ely 
( =1. G- PSTAePSI Ade Rd 


At =CABSEo) 
IF (AD LE. 4, @€-6 > GO TO 243 


WRITE (6. e@@&> 
FORMAT «°@ GCENCMINATCOR IS ZERO . PSIA TS ToC 
10 TO aso 


* BACKSCATTEREC FIELD + a 


BSC = (STi + SED + SSSd*# C5 


LARGE > 


am] 


MANMNANAOANM 


CKL = C4. G@* PIC CKK*eXKD 
Cke = SGRTCCK13+*@. 2546 
ee EBSC = S@RTCSIGMAS WITH PHASE REFERENCED TO FRONT #* 
* EDGE OF WIRE + 
EGSC =CBSC#* CkKe) 
ETTR¢CIFWS= REALCEBSC) 
BRCTFWS = ETTRCIFWS 
ETTI¢CIFW3=-AIMAGCEBSC) 
NICIFW) = ETTICIFW) 
&@@ CONTINUE 
IF «KP WE. 422 GO TO 966 
Pa Fer L = NMIN, NMAS 
SIGMACL® = 14 @ * ALOGLG@CETTRCLI*ETTR CL? + ETTICLI#ETTICL 3 > 
frre CONTINUE 
; WRITECG, $@G@63 CFROCJ), SIGMACT), J = NMIN, MMAR 
SQ@G@ FORMAT ¢ “1FREGQUENCY RESPONSE OF A THIN WIRE *. “¢, FREGLUENCY 
1 Cross SECTION ep eee. CEs t, So eae 30 
968 CONTINUE 
999 CONTINUE 
FETURN 
ENE 
SUBROUTINE EXPI ¢ ARG2, EIRR, EIKI, KPO 
* THIS SUBROUTINE COMPLITES THE REAL AND IMAGINARY FARTS GF THE 
+ EXPONENTIAL INTEGRAL ECs WHERE 
* Eck) = INTEGRAL FROM © TO INFINITY CF EXPCITIAT#OT 
* REFERENCE - FUNCTIONS OF MATHEMATICAL FHYSICS EY MAGNUS AND 
* CBERHETTINGER. FR. or-oe 
KP NE. @ PRINT GUT REAL ANG IMAGINARY FARTS OF ESPONENTIAL 
INTEGRAL 
COMPLEX AIX, ANTGS 166), EI#, FA, F 
FEAL AINTRC16G>, AINTIC1@@) 
DEL = @ GhGee4 
GAMMA = @& SPreis66 
IF ¢ ARG2 .LE. 1. 68-6) Ga Th So 
ABZ = ABS*CARGZ? 
po & J=4,106 
AINTRECT>= & & 
AINTICI3= @.@ 
ANTG¢I> = CMELECG. &, &. &) 
§ CONTINUE 
LE AGE VGE. 25? GU Log 
THIS SERIES USED FOR G@. LT. ABSCKI.LT 415 
** SERIES EXPANSION INVOLVING CICK) AND SICK) 
Ecz>) = GAM + LNCZ) +5UMN CCC -LIeeN eZee SEND OS CC SaeN DEF ACT OSeN DD 
+ JHE HSUM CC CHL I HEN HZHECSHENEL OCC LENE LHP ACT OSH NEL - 
PIO2 2 
INC =1 
FAC = -CARGZ*#29/C2. &) 
AINTR(1) = FAC * @.& 
AINTI¢C1) = FAC#CARGZ/S. @) + ARGZ 
1@ INC = INC +14 


v2 =(2* INC) L-S% 
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1606 


46 
6G 


114 


Ko >1.@ 


AINTRE INC? = FACS + AINTRCINC-1) 
AINTICING) = FACHCARGEXEP159 +FIINTicINC-1) 


AINTRAINCS » 


AER = ABS« 

HERMi= ABS. AINTRCINC-1)) 
HE, = ABSC AINTICTWCd > 
HBLiiL= ABS< AINTI.ING-19) 


IF ¢ ABSC ABR - nERNIe GE Pei wer 7 
IF ¢« AEBS ABI - ABIMNG? . LE. DELS oe Tad 4a 
CONTINUE 


RE C FNS LT. Lee. > GO TG fe 

WRITE 6, 2@G6% 

FORMATS “@SERIES CIC NOT CONVERGE 

NI = INC 

Gad Ta é&& 

NI= INC 

CONTINUE 

EIXF = AINTRE NI? + ALOGCARGO) + GAMMA 
EIXI = AINTICNID - 1. S*O79ES 

GO To 75 


** ASYMPTOTIC SERIES EXPANSION FOR INT FROM INF TO 4 OF 


CEXPC-ITICTI DT 


CONTINUE 

Te « AGS . GE sf? Gi To 39 

SERIES EXPANSION USEC WHEN 18. LE. ABSx> LT. 1568 

ECK) = EXPCIAIECISIN + LCI XD eee + BFACT CIR) eee t 


AIX = CMPLEC@. @, ARGS 3 
FA = 1. G6-°AIX 
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F = FAFA 

ANTGCe2) = F + ANTOCL) 
INC = @ 


ING = INC + 1 

KF = INC-1 

Fo = FeeFePA 

ANTGCING) = F + ANTGCINC-1) 

AB = CABS CANTGCINC)) 

AML =CABS CANTGCINC-1)9 

IF € BESKAB - REAL) -LT. GELY GO TOLLS 
ie « ENG . LE. £9) GO TOLie 


WRITE ¢ 6, 1001 34 
FORMAT ¢ © SERIES DIC NOT CONVERGE “> 


NN = 15 
GO To 74 
NN = INC 
CONTINUE 
EIX = ANTGC NN D*CEXP CRIN) 
EIXR = REALCEIN) u-S$4 
eile MEE ES 
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EIXI = AIMAGCEITS®2 | 
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FUNCTION 
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CONTINJE 


SET TO ZERO IF ARGUMENT ABSee>. GE. 156 
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GO TO 96 
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FORMAT « 
CONTINUE 
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CONTINUE 


RETLRN 


END L-eas 


OR IS ZERO® 


PUREAL FAKT OF EIX)= "915 e5e//9" IMAGINARY PAST OF E(X) 


) 
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L.5  FRUSTRUM - CYLINDER - FRUSTRUM 


The far-field scattering from a frustrum-cylinder- 
frustrum target configuration shown in Figure L.5-1 has been 
formulated using the Ruck-Ufimtsev technique (Ref. 6). 


3 
The expressions of the target frequency response are 
the following: 


v= + ¥mg(1) + 8(2) + 8(3) + (4) + B(5) + 8(6) + 
H 

g(7) + 8(8)} 
where (1) the g(m) are the sum of the scattering due to 


the uniform and non-uniform current associated 
with edge m 


e(Q) 


(2) The upper and lower signs correspond to vertical 
and horizontal polarization, respectively, 


(3) e tC harmonic time variation is assumed, and 


(4) the scattering geometry is presented in Figure 


5-1 
For 0<@<7/2) 
P a % + 
g(1) = aye Bors C(n,) ¥ B(a,,"/248) + 0.5 tan(a, ar, | 
nia, yey 


g(2) = age 2 {4s,,_ [C(ny) ¥ B(ny,4,+0) F 0.5 tan(a,+9)F, | 


= S(a,) JJ, 


= ip re 
g(3) = ae 3 {3ty9_ [ c(ng) F B(ng,@) + 0.5 tan (a5+0)F5 | 


- C(ng) JJo5 


= Py + em 
g(4) = a,o P& cee [ C(n,) + B(ny,@)+0) ¥ 0.5 tan(ay+e)F, ] 


- C(ny) Io, 
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Fig. L.5-1 


GEOMETRY OF FRUSTRUM*CYLINDER=FRUSTRUM 
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s(5) = aye L{ gs, 14 [0m )@, F B(ny, 7/2-0)Qs 


+ 0.5 tan(a,-6)F5Q_ | - C(my) JJp195| 


: ‘i ip + es 3 
g(6) = a,e"P2 {JJ,,, [ c(n,) * B(n,,4,-) ¥ 0.5 tan(a,-0)F, | 


= (ny) ISyy } 
e(7) = + age P3 [35,54 (0.5) tan (27-0) [ 7-1] Qe 
of = ip we z big 3 
g(8) = aje'Pa fas, [C(n,) ¥ B(n,,a,-6) F 0.5 tan(a,-0)F, | 


- C(n,) IIy4} Qe 


For 7/2 @ 


g(1) = aye "1 {J5,,_ [c(n,) ¥ B(n,,7-0-a,) + 0.5 tan(a,+0)F, | 
- €(m) J3p)} Q 


g(2) = a,e 2 tad [ C(ny) Q + B(n5,7-8)Q, + 0.5 tan(a +0) FQ, | 


+ 0.5 tan8(Q,-Q3) - C(ny) 55.99% | 
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AER RBS INLD, EMMY eset sit 


ip3 ha. 
g(3) = age” { Jy9_ [cing) ¥ B(ng,7-a)-6) + 0.5 tan (a,+6)F | 
: C(n3)JJ99 } Q3 
/ EP4 aye 3a — 
2(4) = aye ids. [ c(a,) + B(ny,, 2-9) + 0.5 tan (a+ )F,Q3 | 
- C(n, )IS5y, } 
ip, - 
g{S) = a,e ane [ ccn,) B(n, 6-7/2) | - C(n, IIo, 
(5) = g(6) = g(7) = 0 


where the upper and lower signs in the previous expressions 
correspond to vertical and borizontal polarization, res- 
pectively, and 


sin = 
c(n) = 4 
: n cos 7/n - l 
sin "/n l 
nA) ee se ee res 
‘ n cos ™/n - cos a 


SJimt * [Jp (2ikagsin ) +i Jy (2ka,sin ) | 


= (Nes : : 
JJ5.. [ J, (2ka,sin ) + J, (2ka sin )| 
a 
my * ise — 
a 
a 
n3 =l- 2 
L-73 
. wearer = - = wovcnde sitet " at e re" ee og a 
we 


2k (hy + hy) cos 60 


= -2k hy cos 6 


2k hg cos@ 
2k (hy a hy) cos 6 


Q(2ka, (7/2-6)) 


= Q(2ka (7- 


:) @ 


1 1 


= 2ka 


mes, ia 


- a,-a 
tan L 2 1 
hy 
~-) a,-a 
tan 4 2 
hg 


re -0@)) 


-6)) 


cos (4, +6) 
9\c8¢ 9cos (2, +9) 
4 

5) = pag 0s (41-4) 


cos (a,-@) 


2 


nl 4°49 


as a tan SS 
+2 Tm 
ee oo eee is iit” ge 
fe) 


2n/ = wave number. 


nx 
it 


The preceeding equation can be used in computing the 
first-order scattering from the target; however, it does not 
include the effects of multiple reflection or diffraction. 
The cylindrical surface between edges 2 and 3 is partially 
shadowed for the case of aspect angles between 150 degrees 
and 170 degrees, but the magnitude of this surface reflec- 
tion is small; therefore, this return is formulated using 
the physical boundaries of the cylindrical surface and the 
screening functions rather than the illuminated portion of 
the surface, The screening functions were also used in 
describing the effects of shadowing upon the returns from 
the target edges. 


L.5el Inputs 
The subroutine inputs are read from cards or passed in 
a common block. The parameters passed in the common block 
include: 
NMIN = minimum frequency number 
NMAX = maximum frequency number 
DF = frequency increment (in MHz) 


FC = carrier frequency (in GHz) 


The card inputs are the following: 


MATH 
SYMBOL VARIABLE DEFINITION COLUMNS 
Card 1 ) Aspect Azimuth angle (degrees) 1-7 
ay Al Smallest frustrum 8-14 
radius (inches) 
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MATH 
SYMBOL VARIABLE DEFINITION COLUMNS 
ay A2 Radius of cylinder 15-21 
(inches) 
ay A4 Largest frustrum 22-28 
radius (inches) 
hy H1 Length of first 29-35 
frustrum (inches) 
hy H2 Half-length of cylinder 36-42 
(inches) 
hy H3 Length of second 43-49 


frustrum (inches) 


L.5.2 Outputs 


The data base output consists of four linear arrays, 
EVVR, EVVI, EHHR, EHHI, which contain the real and imaginary 
parts of the vertically and horizontally polarized back- 
scattered fields (in meters) at frequency increments spaced 
DF MHz from NMIN*DF to NMAX*DF, 


L.5.3 Restrictions 


L.5.3.1 Physical Dimensions 


All target dimensions should be large with respect to 
the wavelength of the illuminating field. In addition, the 
target dimensions should be selected such that the angles 

1 and 2, of the frustra retain the basic target shape, 
i.e. the target does not degenerate into a cylinder. A 
restriction on these angles is the following: 


15°<a@4 < 60° 


20°<a, < 60°, 


L.5.3.2 Output 


The output arrays are passed in the argument list and 
a value is computed and stored only in locations NMIN to 
NMAX, 
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L.5.3.3 Azimuth 


The azimuth angle is restricted to the region between 
0 and 180 degrees. In addition the specular azimuths of 0, 
90, 180, (90 -@,), and (90 - &>) degrees should not be 
used. In order to compute the response at these angles, an 
angular offset of approximately 0.01 degrees should be used. 


L.5.4 Definitions of Selected 
Terms Used in Subroutines 


sin = 
1 


————_—_ 


SOREL = CY = re et 


where n= ny) = 1.5 - = 


COEF11 = cos 7/n - cos 


where y = a For @ <7/2 
V _ 

COEF] q = C(n,) + B(n,,7/2+8) For @ <2/2 
Vv _ 

COEF] y = C(n,) + B(ny,7-0-2)) For @ <7/2 


PHASE 1 = 9!P] 


where P) = ~2k (hy + hy) cos 0 


FF ¥V 02 = a,e1P2 (5342. [ cing) + B(nys2,+9)]} ror @ <n/2 


VV 


i ip 
FF yy 09 = ¥ 0,5 tan(a,+®)Fy ae 7 {3,2. | For @ <7/2 


i 
FFW17 = - C(n,) JJ,, age P2 “Far eS af2 


2 


= 7 
TERMIP =JJ)mF = [ 5, (2ka,sin ) +i J, (2ka,sin )| For @ < 2/2 


where m = 1 and the + (lower) sign is used. 
L.5.5 Subroutines Used 
Subfunctions: 
1. Q(X) computes the exponential smoothing function 
2. F(TAUS) computes the special F function 
Subroutines: 
BESL(SI1, XJO, XJ1L, XJ2) computes and returns 
Jo(SI1) in xJO 
J, (SIL) in XJ1 


J,(SIL) in XJ2 
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46. LO 


——- 
am 


SUBROUTINE TARGET 
rr eee 


CEVYR, 


EVVWI, 


EHHR, EHH 


I} 


ie be ST-TA -- FRUSTRA-CYLINGER-FRUSTRA CUFIMTSEY) 
c ———————————— 
COMMCN MOVER, M. NMIN, NMAN, DEL FC. FM, T@ 
is NMIN = MINIMUM FREQUENCS SAMPLE 
' NMAN = MASIMUM FREGUENCY SAMPLE 
C DE = FREQUENCY INCREMENT IN MHZ 
c FC = CARRIER FREQUENCY IN GHZ 
COMPLES TERMIP, TERMIM, TERMZP, TERMZM, TERM4E, TERM4M, 
1 FHASES. PHASES, FHASE4, FFYVOL, FFYVO2, FFYVOR, FFYYO4, 
2 FEWVGE. FEYYaT, FEVVGS, FEYVaS, FRVIAG, FRvYad, FRYYL2, 
= FEVYL4, FEVVLS, FEWYLG, FEVVL?. FRVOLS, FEYVLS, FRVY2O, 
4 FFVYVeo, FFHHG@1, FRHHYS FFHHG?, FRFHHG4. FRFHHGS, FRRHG@E., 
5 FFHHGS, FFHHGS, FFHH1G, FFHH14, FFHH12, FFHHiz, FFHH44 
& FFHHLG. FFHH17, FEHH1&, FFHHLS, FFHH2G, FFHH24, FFHHSe. 
7 FFHH, CEYY, CFHH,  FTAUL, FTAUZ, FTAUZ, FTAUA, 
& FTAUE. FTAUS FTALS 
COMPLEN F 
L 
DIMENSION EVYRESi2>, EVYICS12), EHHRCS129, EMHICS412) 
ind L 
: C + + * ALL DIMENSIONS ARE IN INCHES AND FILL ANGLES ARE IN 
C 
READCS, 4000) ASPECT, A4, AZ fd, Ht) H2, HE H4, HS M4 
1000 FORMATCOF?. 2,12) 
WRITE (6, 101@) AL, Ht, AZ, He. Aa, HE 
1010 FORMAT ¢ 21H ° FRUSTRA-CYLINGER-FRUSTRA °) 
1 SON * UFINTSEY-RUCK FORMULATION ¢./ 
2 é ft « -FRode BRO MER OF el ayy 
3 GW «AP = . Fis ON Hee LES. Sc 7 
4 SH Aa = F672, OH HE =|, F803,73 
c 
c = 44. sazesare 
PI = 2 141592e5rsears 
SPI = SGRTCPID 
{ TR * PI ¢ 180.8 
RTD = 480. @/PI 
We = 2.@« PE * FC 
#KO@ = WCC 
WekKG = SKGG+KKOG 
wEKOAL = NOK@RML 
MEKOAZ = ¥2KAKRS 
K2KGA4 = K2EKOKAA ‘ 
C 
THETA = ASPECT * DTR 
STHT = SINCTHETAD 
CTHT = COS¢THETAD 
TANATT = STHT ¢ CTHT 
Cc 
SHADOW = (Ad - AZ) ¢ CHE + HE + HDD 
SHADOW = ATANC SHADOW) 
ALPMAL = ATANCCAZ-AL)/HD) 
ALPHAS = ATANCCA4-AZ)/HED 
¥1D = ALPHAL#RTO mn be | 


a 


* * 


PHASE, 
FFYYAS, 
FFYY413, 
FFVYV21, 
FFHH?, 
FFHHLS, 
FRY, 
FTAUS, 


CEGREES + * 


- 


nl 


<G16 


as 


Keb = ALFPHAS#kTO 

NEC = SHACOMERTO 

WRITE C6. 2010) ASPECT, X10. 820. #20 

FORMAT ¢ 18H@ ASPECT ANGLE =, FE&. 5, -¢; 
4tH ALPHAL = 7FPE E77 14H ALPHAS = 
19H SHADONWCALPHAZ) = Fa. £3 

SAL = SINCALFPHAL? 

SAS = SINC ALPHA] 

ALFPT = ALPHAL+THETA 

AFLMT = ALFPHAL-THETA 

AePT = ALPHAS+THETA 

AéeMT = ALFPHAS-THETA 

TANAPL = TAN ALET> 

TANAPe = TANCASPT 2 

CALPTS = 2. #COS* ALFETO “Sri 

CHGFTS = 2, eCUSCREFT) “She 

MND = 1.5 - (ALPHAL FId 

Me = 41.6 + CALPHAG Fil 

No = 1.8 - (ALPHA? Fl} 

‘Na. = 4.5 + (CALPHRe PI? 

CPONL = COSCFICKNIL? 

CPONE = COSCFISCRNSD 

CFON: = COSCPICRNED 

CEOs = COSCRICKN4) 

TERMG@L = CSINCFI ¢ SNLD) aM! 

TERMGS = CSINCFI ¢ RNBO) HANS 

TERMNG@S = (SINC KNEDD ¢ RNS 

TERMG@4 = CSINCEI &¢ XN4)) eNG 

COEFNS ARE CONK) TERMS 

COEFNL= TERMG@L1(CCRONL - 1. > 

COEFNE= TERMG@Ye“CCRONE - 1.3) 

COEFNZ= TERM@ECCRPONE - 1, > 

COEFNM4= TERMG4“CCRONG - 1.5 

IF CASPECT .&). $663 (Gh Te 26 


CIF FRACTION TERMS ¢ CCND-“+BCN, FHI = 
BG 8 ls Vs 


COMPUTED HERE FOR 


CUERI. = 2.6 ¢. %-CRURS = 
COEFLeE © 4 6 € CFOUNE 

CUEFIS = 1.0.0 € UPONS = 
COEFI4 = 1.6 ¢ <« CPON@ - 
COEF22 = 2.60 7 € CPONL - 
COEF2Z2 = 1.8 /7 € CRONE - 
COEFe4 = 1.6 7 ¢€ CPON4G - 
COEFAY = COEFNA - COEF1I1 
COEPIH = CHEFNL + COEFLi 
COEFeY = COEFN2 - COEFI2 
COEFEH = COEFNe + COEFAe 
CGOEFZY = COEFNS - COEFA4= 
COEFSH = COEFNZ + COEFA= 
COEF4Y = COEFN4 - COEF14 
COEF4H = COEFN4 + COEFI4 
COEFSY = COEFN1L - COEF21 
COEFSH = COEFN1 + COEFEe1 


COEF RX 


,FE 


TERNS 


COS¢¢ THETA+THETATFEI OC RNG 2D 
COSC CALET + ALPT) AHNSD > 
COSCCTHETA + THETA) “XNZ)) 
COStCAZPT © REPTS ANN4)) 


COSCXTHETAFT HET AP ior 
COSCCALMT + AINTo) c 
COSCCASMT + ASMT) e 


TERMG4 
TERMG4 
TERM@e 
TERMGe 
TERNGS 
TERNMGS 
TERNG4 
TERNG4 
TERMNGL 
TERMGL 


PA en eae eee & 


*N1)> 
wNe2))> 


4 
‘ 


sN4)) 


COEFEY = COEFNG - COEF22 * TERMGEe 
COEF6H = COEFN2 + COEF2@e * TERMGe 
COEFFCY = COEFN4 - COEFZ4 *« TERMG4 
COEFTH = COEFN4 + COEFS4 + TERME 
L. 
TANANY = TANCAIMT> 
TANANG = TANCARSMT 3 
CALMNTS = 2. #COSCAHLMTI SAS 
ECOSCAEMTI “SAS 
+ ALMT , 
* CSHADGH - THETA) ! 
c 
C CIFFRACTION TERMS ¢€ CCND-¢+BRON, FHI) = COERFRM TERMS 
C COMPUTEG HERE FOR THETA GT of 
1@ COEFZ4 = 4.6 ¢ CCPONL = GUS¢ 2. Geer DFAT RNS OD 
COEF <2 = @ ¢ CCPRON2 - COS¢e. Ge FI -THETAI KN) 4 
COEF ERE = 4.6 ¢ CCPRONS - COSCS, GeC(PI<AZPT) CRNS)) 
COEFR4 = 4.6 4 CCPON4 = COSKC2. G# THETA - ©. GePID CKN4G DD 
COEF 44 = 7. @ ¢ CCRONG - COS¢( 2, GHTHETA = FI 17RN4DD 
C 
COEFLY = COERPNE —- CHEFSL + TERMEL 
COEFLM = CERNE * COEFSI * TERMEL 
COEFeV = CUEFNe - COEFSe * TERMG@2 
‘ CUEF ZH = CUBRPNe + COEF ¢ TERMS 
COEFZY = COEFNS - COEF * TERME 
COEFSH = COEFNE + COEFSS * TERMG@S: 
COEF4Y = COEFN4 - COEFZS4 * TERMG@d 
COEFGH = COEFN4 + COEFE4 + TéERMad 
COEFSY = COEFNG - COEF44 * TERMG&4 
COEFSH = COEFN4 + COEF44d + TERMAd 
\ 
GOL=00 82K HAL # 0 PI-ALFETS ) 
Q2=OcnekKGAs * CFI - SHADCGW - THETA»: 
GSC KheKG@AS * CFPI-AZET) ) 
26 CONTINUE 
C 
Do 9@6 I = NMIN, NMA 
’ 1 SD acl 
MW = (2. .% Pl. KE & GFX #2 PEGGIG 
aKa IN ae eine 
BRE G2AC = KK * CCTHT + CTHT) 
{ 
TAU = #K@ + AL * CAAPTS 
TALE = \K@ * Ae * CALFTS 
TALIS = KKGO * Ae « CAZPTS 
TALI4 = KKG * Ad * CAZFETS 
FTAUL = F e TALIA 
FTAUZ = F (TAU) 
FrAUs =. F CTALIZ) 
FTAUd = F CTRALI4D 
C 
SES = 2.@ * KK@ * AL * STHT 
Sl2 = 2.0 " KK@ * A@ * STHT 
S14 = 2.60 * KK@ * AG * STHT 


CALL BESL (S11 KJGH1, KJIe4, KJ2e4) L-8/ 
CALL BESL (SIe@, KJGK2, KJ1Ike, KJeXxed 
CALL BESL (S14, wJ@K4, KTITKA, KJeK4) 


TERMNAF = CMPLECKIGKL, HILK19 
TERMIN = CONJG*¢TERMIF) 
TERMNEF = CMFPLYCN IONS, BILK SD 
TERMEM = CONJGCTERMEF) 
TERM4F = CMPLE CN TOR4, MIIK 4D 
| TERMd¢M = CONTIGS TERM4FS 
TERMS = KJ@X4i + RUSH 
TERME = XJ@Ne + | re 
TERM? = 8JGK4 + RISK 
PHI4 = NKG@2C * CHItHe? 
FHI e = “NKG@eC * CHE?D 
FHI = FHIe 
FHI4 = eEGEC * CHEtHD3 
PHASEL = CMPLEACCOSCPHIL), -STNCPHIED > 
PHASES = CMPLACCOSCPHI2), -SINCFHI2) ) 
PHASES @ CHPLACCOSCPRTE 2S), SCRE yD 
FHASE4 = CMPLRCCOSCPHI4?. SINCPHI 4S 5 


LF SHSPECT “690 SGeo re aie Tih ys 


TAUS = Ka * AL + CRAMTS 
Tae = KG * Ae * CALMTS 
TAU? = a Sr & CHLMTS 
Pause 2 Ke * Ad *« CAZMTS 
FTAUS =< oF THLIS 3 

ae FTALE = fF CTRUE 3 
FTE = F CTALP 
ETS oe TRUE) 
FEY Vil = AL * TERMIM * COEFIY + PHASES 
FF HeRaL = Fd * TERMIM * COEFLH + FPHASEL 
FR u's = Ae * TERMEM + CUEFGY + PHASES 
FFHHie = Ae * TERMEM + COEF@H * FHASE? 
FFYVGE = Ae * TERMeM * COEFSY & FHASES 
FFHHG@2 = 42 * TERMEM * COEFSH + PHASES 
Frvved = Ad » TERM4M * COEFGY «# FHASES 
FFHHEd = Ad * TERM4SM + COEFS4H + PHASES 
FFVWES = At * TERMIP »« COEFSY * PHASES 
FFHH@S = Al * TERMIP * COEFSH * FHASEL 
FFVV@6 = Ac * TERMEP & CORFEY « FHASES #* Ee 

- FFHH@E = Ac + TERM2P * COEFEH + FHASES * GE 

FFYVYGr = #4 & TERM4P + CIERPTY + PHASES «# Ge 
FFRHG? = H4 + TERMGP * COEFTH + FHASE4 * GE 
FFWVG@E Al * TERMIM + TANAPL « « @& SD * FTALL * PHASED 


FFHHES 
FFYVaS 
FFHHaS 


~FRYWY EE 
Ae »* TERM2ZM 
~FRFYYVas 


= 


TANFAIPA # ¢-@ 5) # FTAUE * PHASE? 


FFYWW16 Ae + TERMEM * TANAPe * ¢ @ SS) * FTAUS * PHASES 
FFHH16 -~FFYYV16 is 

FRFVY11 Ad * TERM4M * TANAPS * ©-@. 5) # FTAUd * PHASES 
FFHH1i1 =PPV VALS 

FFWW1e Al * TERMIF * TANAM2 * ¢ @.&) * FTAUS * PHASEL * G6 
FFHHLe se al oe ih 8 

FRYYV4< Az * TERMZP # TANAML # ¢-@& 5> * FTAUG * PHASES * G6 
FFHH12 ~FFVV42 


FFWY14 
FFHH14 


Ae * TERMZF 


+ 


TANAMeG #¢€-. S)#C4. -FTAUZ)*® PHASES * G& 


SF SB He FEB NR UD HR 


~FFYY14 L-92 


28 


FFYY4S 
FFHH1S 
FRFVWOIE 
FFHH16 
FRA? 
FFHH1i?; 
FFYYIE 
FFHH1IE 
FRFYVY49 
FFHHLS 


now ou 


FRVWOO 
FFHHZG 
FRYYoL 
FFHH24 
FRYYoo 
FFHH22 


Hs me 


vou on 


FF YY 


LP Ol 


FFHH 


tad Po 


oo Td 46 


CONTINUE 


FEY Yad 
FFHH SSL 
FRWY ie 
FFHHHe 
oe 
FFHHES 
FFW Yad 
FFHHE4 
FRYVES 
FFHHES 
FFWY EE 
FFHH@E 
FFYVA? 
FFHHE? 
FFYVEE 
FFHHGS 
FRFYYW@s 
FFHHaS 
FFVY414 
FFHH11 
FFVY12 
FFHH12 
FRVV43 
FFHH1: 
FFYW14 
FFHH14 
FFYY45 
FFHH1S 
FFVY16 


ae ewe SE EES EH ES EBA He eR Re eee 


Ad + TERMGP & TANAMNG + ¢-@ 5) # FTAUS * PHASES * 6& 
-FRWW4A5 

-Al * COEFNL * TERMS © PHASE. 

FFYWIE 

-Aic * COEFNG *# TERME # PHASES 

FRY? 

-Ae * COEFNE * TERMGé * FHASES 

FRYY SS 

-H4d + COEFNd + TERM? + PHASES 

FRYW AS 

-Al * COEFNL * TERMS + PHASES 

FEW oe 

-A@ * COEFNG * TERME + PHASES * G6 

FRY Wod 

-Ad * CTOEFNS * TERM? + PHASES * OF 

FE Wwe 

FRWUWEI + FRYYEo + FRYWOE + FRVYGd + FRY YOS + FRY VGE + 
FRYWO? + FEFYWaS + FRY + FRY WAG FFYY41 + FRFYY42 + 
FRWWAE + FRYY14d + FRYY1S + FRYYLE + FRYVAT + FRYVYLS + 
FRVW4Ao + FRVV2G + FRYYod + FRYWoe 

FFHHG@t + FFHHGS + FEHHGE + FRHHG4 + SPHHES + FFHH@E + 
FREHHE? + FREHHeS + FREHHOS + FRHHIG@ + FREHH12 + FEHHL2 + 
FRHHL= + FEFHHid + FRHHIS + FRHHIE + SP YHA + FRHHLS + 
FFHHIS + FREHHOG + FREHHEd + FRHHeS 

Ad #* TERMIiM * COEFLY # PHASE. * G4 

fil + TERMIM + COEFLH + FHASEL ® 34 

fhe * TERMEM * COERSY + FHASES # Ge 

Ae * TERMSM * COEFGH * PHASES * Ge 

AZ * TERMZM * COERFZY # PHASES # G3 

A2 * TERMeM * COEFZH * FHASES »* GF 

Ad * TERM4M * COES4$Y # PHASES 

Ag + TERM4M # COEFGH + FHASE4 

Ad * TERMN4F * EOERSY * PHASES 

Ad + TERM4F + COEFSH # FPHASE4 

Al * TERMIiM # TANARL ¢ © & S) + FTAUL * FHASEL * O41 
~FFYVOE 

AS * TERMGM * TRNAPD * ©-@ 59 # FTAUS # FHASES * OL 
-FRYYE? 

Ae * TERMZM * TANAPS * ¢ & 89 © FTAUE + PHASE: * @3 
-FRFY YE 

Az * TERM2M * TANATT * ¢ @ S) * PHASES #*¢(@2-Q3) 
-FR UVES 

Ad + TERM4M * TANAPS * © -@ 5) * FTALIG + PHASES * GF 
-FRYV4O4 

-A1 * COEFNY * TERMS + FPHASEL * G4 

FFYW42 
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14006 FUNCTION C(Z) RCST 424 
1401 c G(2) = O.58(1 + ERF(Z)) RCST 425 
1402 c * ERF(Z) IS EVALUATED USING A RATIONAL PCLYNOMIAL APPROXIMATION RCST 426 
1403 c * REFERENCE (HANDEK MATH FUNCT RY ARRAMOWITZ AND STEGUN, RCST 427 
1404 c * SECTION 7.1.26) RCS7 428 
1405 c RCS? 429 
1406 IF ( Z.GT. 2.) GO TO 10 RCS7 430 
1407 IF ( Z.LT.-2.) GO TC 20 RCS7 431 
1408 AZ = ABS(2) RCS7 432 
14049 P= LeO/(1.0 + «47TO4T*AT) RCST 433 
1410 Y = 1.0 — P*®( 3480242 -— P*(.O9587T9E — .747E55H*P) DEEXP(-AZ#AZ) RCST 434 
1411 IF (2) 29496 RCST 435 
1412 2Q = (1.0 - Y)/2. RCST 436 
1413 RETURN RCS7 437 
1414 4Q= .5 RCST 438 
1415 RETURN RCST7T 439 
1416 6 Q = (1.0 + Y)/2. RCST 440 
1417 RETURN RCST 441 
1418 10 0 le RCST 442 
1419 RETURN RCST 443 
1420 20 Q@ = 0. RCST 444 
1421 RETURN RCS7T 445 


1422 END U- Alen RCST 446 


L.6 CYLINDER-FRUSTRUM COMBINATION 


The far-field scattering from a cylinder-frustrum com- 
bination target shown in Figure L.6.1 has been formulated 
using the Ruck-Ufimtsev technique (Ref. 7). The solution 
can be represented in the following form: 


e(@) = * Va {g(11) + g(12) + g(1) + g(2) 


+ g(3) + g(4) + g(9) + g(5) + g(6) 
+ 8(7) + 8(8)| 


where g(m) represents the sum of the fringe wave scattering 
and the physical optics response associated with 
edge m and the upper and lower signs correspond 
to vertical and horizontal polarization, 
respectively. 


With the use of the diffraction coefficients at the 


concave edges (F #1), the expressions for the g(m) are the 
following: 


For 0<90@< 7/2 


+1 


ip 
g(11) = aye P21 / 53,4, [c(n,) F B(n,, F - @)]- c(n,)IIyy |} 


+1 


i 
g(12) = aye Pll {53,,_ [c(n) ¥ B(n,, F + 0) - C(n_)IS,, } 


_ip - = 
g(2) = ase 2 | JJy9_ [C(ny) ¥ B(ny,0+a) ¥ 0.5 tan (a+0)Fo] 


- G(n2)IJp2 | 


Fig. L.6-1 GEOMETRY OF CYLINDER-FRUSTRUM COMBINATION 
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ip ~ = 
g(4) = aye te | (ny) + B(n, , +a) + 0.5 tan (a+0)F,| 


. (ny) II, | 
= a.erP9 ? 
g(9) = ae P9{53,,. [c(n,) * B(n,@)] - C(n,) 335, | 
g(1) = a,e'Pl JJ}1- [c(n,) + B(n, ,9) + 0.5 tan (a+0)F; | 


g(3) = aje'P3 {33,5_ [c(ng) ¥ B(ng,0) + 0.5 tan (a+0)F5] 


i C(n4) JJ, | 


¥ 0.5 tan (a-9) [1 - F, a, ] 


a ip 
g(5) aje 1 Meat 


g(6) = aye P2{ Wiss [ccn,) + B(n,,«-0) + 0.5 tan (a-0)F ¢] 
- C(ng)JJ99 } 26 

g(7) = aye P3 J5p4{ ¥ 0.5 tan (a-0) [1 - F, ]}9, 

g(8) =a eP4 | IS 4 [ c(n,) + B(n, ,a-0) + 0.5 tan (a-0)F | 


4 


- C(n,)II,, 1A 


For 1/2<O<7 
g(5) = g(6) = g(7) = g(8) = 0 . 


g(9) = a,e°P9| JJ4,.| c(n.) + B(n, ,9)] - c(n_)5J5,} 


g(10) = a,e'P9 | 2 [e(n,) 7 B(n, 6-3) | ‘ C(n,) IIo4| 


g(4) = a,e Pe { J5y4. [c(n,) * B(ny,7-) ¥ 0.5 tan (+0)F,Q, | 
- C(n,)JJ,, | 
g(2) = age? | JJq9- [e(ng)Q + B(ny ,7-8)Qo 


F 0.5 tan (@+0)F2Q | - Cy) IT9 9% | 
g(1l) = aye PLL (53, ). [c(n,) + B(n,,7-0) | - C(n,)JJo} Quy 
g(1) = aje'Pl {55,,_[C(n,) ¥ B(ny ,w-a-0) 
+ 9.5 tan (a+0)F, |Q, 
-C (ny) 35p4% + 0.5 JJ,,_ tan Q,, [1-e,]} 
g(3) = ae > { 33yp_ [C(ng) ¥ B(n3, 7-2-0) 
+ 0.5 tan (a+0)F,]Q, 


- C(n3)II4%3 = 0.9: Sy Sane 8, [1-2, ]} 
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+ ee eee paren ee es 


ee Sek: 
> ey OR 


where the upper and lower signs in the previous expressions 
correspond to vertical and horizontal polarizations, 
respectively, and 


sin = 
C(n) = = n | 


cos Z -l 
n 


sin = : 
B(n,¥) = 7p aa? 


cos =~ - cos 
a 


ISimy = [Jo (2ka,sine) + i J, (2ka,sind)] 


a . + : 
IIo. [ 5, (2ka_sine) J, (2ka_sin@)] 


7. * 3/2 


Pi, * -2k (hy + hy + hs) cos@ 
Pp, = -2k(h, + hy) cos® 

Po = -2k hcoug 

P3 = 2k hy)cos® 

bad 2k(h, +h,) cos@ 


p5 = 2k(h2 + h3 + hy) cos® 
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"e = Q(2ka e) v4 -0)) 


Q = Q(2ka (tT-a ‘ - @)) 


r? "hts 2ka csca@ cos (a+) 
2 
3 
4 


FNNeE 


7? > 2ka 1\ csca cos (a-@) 
6 2 
7 2 
\8 4 
-a 


a= tan"! 21 = tan ——< 
1 hg 


The expression of the g(m) modified for use in the 
second formulation (F #2) of the scattering from the 
concave edges (1, 3, 5, 7) are the following: 
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a 


For @< 7/2 


g(1) = + aye Pl gy,,_ (0.5) [tan (@+0) (F,-1) + tan 9] 

g(3) = + ase 3 J5zp_ (0.5) [tan (@+0) (F3-1) + tan 0] 
For 0 >7/2 

ci) <4 4 ae [ tan (a+0) (Fy-1) + tan 0] Qy 

g(3) = + ase P3 3J,>_ (0.5) [ tan (a+0) (F3-1) + tan 0] Q). 


The formulations of the basic scattering from Target ST-2 
involved only first-order diffraction and were obtained by 
using the Ruck-Ufimtsev technique. Such higher order scatter- 
ing as multiple diffraction among the edges, multiple reflec- 
tion between the surfaces, or scattering mechanisms involving 
both a creep path along the surface and edge diffraction 

have not been modeled. However, by comparing the signatures 
computed from the basic formulation with wideband measurements 
data, the significance of the higher-order scattering mechan- 
isms could be determined and a decision can be made concern- 
ing which is the better formulation of the basic scattering 
from the concaves edges. 


L.6.1 Inputs 
The subroutine input parameters are read from cards or 
passed in a common block. The parameters passed in the 
common block include: 
NMIN = minimum frequency sample number 


NMAX = maximum frequency sample number 


DF = frequency increment (in MHz) 


FC = carrier frequency (in GHz) 


The card inputs are the following: 


MATH 
| SYMBOL VARIABLE DEFINITION COLUMNS 
Card 1 8 Aspect Azimuth Angle 1-7 

(Degrees) 

ay Al Radius of smallest 8-14 
cylinder (Inches) 

a, A2 Radius of middle 15-21 , 
cylinder (Inches) 

ay, A4 Radius of largest 22-28 
cylinder (Inches) 

hy H1L Length of smallest 29-35 
frustrum (Inches) 

hy H2L Half-length of middle 36-42 
cylinder (Inches) 

hy H3L Length of largest 43-49 
frustrum (Inches) 

hy, H4L Length of largest 50-56 
cylinder (Inches) 

he H5L Length of smallest 57-63 


cylinder (Inches) 


L.6.2 Outputs 


The data base output consists of four linear arrays, 
EVVR, EVVI, EHHR, EHHI, which contain the real and imaginary 
parts of the vertically and horizontally polarized back- 
scattered fields (in meters) at frequency increments of DF 
MHz from NMIN*DF to NMAX*DF, 


L.6.3 Restrictions 


L.6.3.1 Physical Dimensions 


All target dimensions should be large with respect to 
the wavelength of the illuminating field. The target is 
further restricted such that two frustra form the same angle 
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-l a5-a 
with, the cylindrical surfaces, i. e., @ = tan A 
= = ‘ 
tan er ae 
h3 


In addition, the basic shape should not be distorted by 
choosing @ too close to 90 or 0 degrees. A value of 30<a<80 
should be maintained. 


In determining the shadowing of the various surfaces 
and edges the target geometry was further restricted such 
that %3 < 4 where these angles are defined by the equa- 
tion. This restriction applies only to the use of the 
Q functions in determining the amplitude weighting of the 
scattered field terms due to shadowing. 


L.6.3.2 Output 


The output arrays are passed in the argument list and 
a value is computed and stored only in locations NMIN to 
NMAX . 


L.6.3.3 Restrictions 


The azimuth angle is restricted to the region between 
O and 180 degrees. In addition, the specular azimuths 0, 
90, 180, (90 -@) should not be used. In order to compute 
the response at these angles, an angular offset of approxi- 
mately 0.01 degrees should be used. 


L.6.4 Definition of Selected Terms 
Used in Subroutine 


Dl = cos @/n 


Dll = sin 77/n 
ar 1 
- cos Wn-1 
where n = a is ns = 1 -& 
{y}l = C(ny) ¥ BCn,,@) for @ <7/2 


PS1 = JS44_ = JSy nF = [ 49 2ka,,si 6) + i J) (2ka_sin@) 


where m = 1 and the upper (-) sign is utilized 


iL = J 
CCl C(n,) IS, 


CTl = 0.5 tan (@+ @) Fy 
V = 
ct yf! = ajetPy| JJ55. [een + B(n, ,8) + 0.5 tan (a+0)r,| 


: C(ny)II54 | for @ < 7/2 


3 = [3 (2ka_sin®) + J,(2ka sin)] 
fe) m m 


L.6.5 Subroutines Used 


Subfunctions: 


Q (X) returns the value of the exponential smoothing 
function 


Subroutines: 

1. BESL (XCx, BxO, Bxl, Bx2) computes Bessel functions 
of first three orders for real argument Xcx and 
returns 

Jo (XCx) in BxO 
J, (XCx) in Bxl 
Jy (XCx) in Bx2 
2. DIFFC (VX,HX,NX,DX,DX1,EX,PHI) computes 


tee 


i C(N) + BCNX, PHT) 


The inputs are NX, DX, DXl1, EX, and PHI and 
outputs are VX, and HX, 


3, FTG (TSx, FTx) computes the Special F function 
using TSx as argument and returns value in FTx, 
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ad same 


6aaL 
A 
£ 


SUEROQUTINE TAF 
a) RESPONSE 


GET CEYWR, EVI, EHHR, EHH] 3 


OF TARGET 


ST-2 


+e 


++ COMPUTED UTITYCING THE RUCK-UFIMTSEY EQUATIONS ## 


COMMON 
NMIWN = 
NMA = 


MOVER, 
MINT MUM 
MAS IMUM 


M, WNMIN, 
FREQUENS 
FREGUENL 


NMAS, ODF, FC, 
Y SAMPLE 
Y SAMPLE 


Fil, TS 


DF = FREG 


pitepmu 
SCM. 7 


INCREMENT IN MH2 


FC = CARRIER FREQUENCY 


IN GHE 


PSi, Fara) PSt PSLP, Fo2r 


COMF S 
C4, C3 Ce. 2245. PTs 
Tis CT2 Ct4, OG. CTS). CF 
COMPLE CW4i, CHI, Ca, CHE, 
6 Ho» CGV?, CH, CVE. CHE 
é Hie CY, CH CVA, CHA 
REAL Ni. Ne, NC. ITM 
TENSION EVVRCL>, EWWVIC1), EHHEC 2), 
FROGRAM CONSTANTS 
PI = S. 142892652569 79 
PIS = FIeP] 
PigGe = FI/'e 
HF = Se PI 
S = SORTCFI 
. = 180. “FI 
(‘TR FI/LE@ 


WRITE (6, 60 
FORMAT: oP af 
rH Ale = 
‘H H4 = 


H12L = HAL 

Hi25L = HAL 
H23L = HEL 

H224l = HeL 
Agel = AZ-AL 
Ade = Fd-FZ 
KL = 


ATANS 


SH 4 PRC 


rH HSFE 


fH Bi 


+ Hel ¢ 
+ Wot 
+ ZL + 


Ael/Hil> 


ASPECT 
GRAS 


INFUT PARAMETERS. 
CT ANGLE = » FE. ¢€ 2 
a4 H1iL, AL, Hel, Ae. He 

= ,FGE.4 ,7H AL = 


rem RS ;FRL4¢ > TH 


»eH HS  ,Fa 4 ) 


HSL 


H4L 


L=\or 


L,A4, H4L, Hee 
»FE. 4.¢ 
Ad = . FS 4,¢ 


Ce, 
Fie, Firs 


CHd,. CVS, CHS: 


4, CHLG, C¥11,CH11, 


EHRI <4) 


»HEL, H4L, HEL, M4 


Atte Hee LP 


oO 


" 


Seere 
n 


— 
4 
<7 

—_ 
te 


10 = 
20 = 
=i 

= & - 
4D = 


& CI 
C& [ 
CIFFC 
CALL 
CALI 
CALL 
rept 


“4 
ATAN( 
BT AN 


Agee 
Aad 


CHEL+HZEL))d 
HIL+HSL > 


= Q¢ 


6, 6865 1D, Med. X20, K40 
T SH Xi © FF. 3 GH NZD = Fr. SE 
fH C30 = , FT. 3 gH ‘> = 2 F?.S 
4. =-Xi7 Fr! 
4. #XL/PI 
4 c 


1 bi-4 
N i 
' * "4 
x N Me 
PI/t NE 
N<TH 
r +Fid 
' Fa 
h PT) *@. 5 
PTO/SINGHS 
4 GT Fla2) Go C 2@ 
tae. 
& S#* TANS KMT 
COSCHATISSINGXSL) 
KeAeeCNe- TH) 
Ri 2A4K&CRE-TH) 


CZ8) 
COMPUTES 

DIFFCC Wie, 
DIFFC ¢ 
DIFFC ¢ 
HIRE ¢ 


CONI-“+E ON, PHI) TERMS RETURNED 
Hi2,NC,0C,0C4,€C, Pide-TH > 
NC, 0C,0C1,EC, FPlQe+TH Dd 

N4, 04,014, €1,TH D L-/o8 

N2,D2,074,.F2. XPT ) 


AS VEL HK 


W41,H414 
V4. H4, 


"2. He, 


eaeembniineehane tt oe 
tes 


- 


tb ee 


ef 


CALL DIFFC ¢ ¥9,H9 .NC, DC, DC1, EC 
» HE 


CALL DIFFC ¢« ¥6 
GQ TO 34 
CONTINUE 
PIMNTH = PI-TH 
e2 = AKSAS*CPIMKE-TH) 
PIMKPT = FI=Ki-TH 
21 = AKGALHPINE PT #S 
e114= AKCAL*(PIMN4-TH) 
@11= @¢211) 

n 


hi Pa Papa 


fe = « 
@i = @ 1 
CALL GIFFC WO, HS,NC. OC, OC41,E 
ALL PIFPE § ‘48, H1iG, WO, O0C, bC4, 
CAL DIFFC Wo. He Ne. Be a Beg = 
AiL OIFFE © ¥VaisHds Nts Ot, C2428 
ALL BDIPFCCVit, Hil, NC. OG. OCT, EC 
ONTINUE 
5 & I MRED NM. ON 
| i | 
AK =AxGek 1 *DF 
AKe = & &+Fik 
ARSAL = AKe*At 
KeAS = eS 
ad AS Air oehd 
; 
Gis $ aS 
Ces aes. F 
4 Sx ST 
ALi BE Ci, 813, 614, 642 
Li we Bot, Bee 
E C4 E4i.bqes 
Fit { BE TERMS a 
PS CMF =B44 
CHF -Es i 
Po4 iP L Bd) 
PSLe c 
PSer = iM 
Poop = , 
PHASE TERM USING LENGTH ALONG Ax 
PCL = ~fike@+HLeElLeTH 
PCS = ~Fke-Hel*eCTH 
PC4 = AKZ*4Z3L*CTH 


POS = AKZeHE TAL ACTH 
2eHi2SL*C7H 


, 
i 
— 

" 
' 
xz 


C1 = CMPLAECCOS<FC1L), 
é = CMPLR‘SCOS¢FC!), 
= CMFPLECCOS*FC4), 
CONTIGS C2) 
CMFLERCCOSCFCS), 


nanan 
rtp & 
— 

now 


CL = AL*SPitCt 


SING POA) 
SINCPC2)) 
SINCPC4)) 


SINCFCS)) 


TH ) 


Ne, be, Bb21,€2, «MT d 


TH-= 
rTMT 
PIMNE 
PTATH 


CNEL ACCOSCPOC119, SINCPC119) 


3 = AS*SPI*€C3 

= AeeSr lee 
4 = AdeSPIecd 

= FidaSPI«c! 
Li = ALYSPI#Cd4 


CAUSTIC CORRECTION TERMS 
CCLi = EC#OCL*(Cb1G+bi2d 
COL = EL*¥DA1* Cb1G+B1 2) 
GE = E2eDedaechoar+bed) 


ZS = EL*Dl eC bo+be2) 
4 = Eee eC bs +b 2) 
g = EC#OCie (he G+ aed 


CC 
es 


= TSPeAK2 
FLeTSPF 
Foe T SPE 


To4= AqdaeTOPe 


robe hy 


7 m4 - 7 
r u Pode Tt 
= Ae ET 4 


a ee - -f 
# pr | me 


te 
Re 

1 
4-4 
- 


ETG C. rT? 
E L4CT So ePSe ~- CCR 
» + 4-CTSIS*PS2 - COED 
f2-CT2I8PSEe = CC2> 
ra Het+CT2o*PSe - CCe2 
AL r 1TS4, FT4) 
T CeerT 4 

Cee eve-CT4i#P Sd - COS 
b tC CHeeCT4#FP Sd ~ CO) 

Coe cV¥oerPsd - COS) 

CHS = CS#eCHSERSd - COO) 


C¥VLLECVLe4+C VL + CVeeCv sel de ly 
CHL4I+C +CHe+CHE+Ch4+CHS 


= 
~ 
&2 
$s 
Lan 
on 
— 


BF 4 26 LE. 2.088) FO Bot 


SMM = TSMAAKE 

§ =TSHM*AL 

Sé =TSMMEAZ 

@ =TSMMeRd L-\\0 


TGCTSE, FTG) 

Té = CKB*XFETE 
Cee CVe-CTEIHFSS 
Cok (CHE FCTE) HP SZ 


= 
ne 
— 
” 


FTG Se ee 
PTS = CKORCL. -FTSIePSiPHCL EGE 
; = =F 
} = +F 
& = CkK@e¢4 - ee Ware 
eTs 
PTA Ts FEE 
= ie HE TE 
2 d+ -CTELAPSAP - Coda +b 
Lo t *U Ta KR Sd P = COd dp whe 
-—-EYS + 3 & * iy 4‘ é 
= + roe + Lf + Fe 
’ cis & 
me ss 
~ ey 
sup 
. e 
4 te 
. se ee ty , 
t inh tr rife Fed 
+ © - 
Ly + f - £099 
1G = CoacVibeP Ste Be 
Le Se CHLOE SEP ~ COS 
= f Shae Pad C4 
+ Serad - Ce 
co + CW4E #044 
i + CHLG +04 
ce LE ~e Go Th Sau 
Yi. ® cia LePSt = CCL etd 
{ = i4 Hi L#Prs - Cli +@42 


Ve * CY GR PS 
Che * C2* CH4eFS 


fe OE 
ve c 


ao 
Pano 
* 
a 
m7 fis 


EV = CY * CVLL €© C¥2 

CH = CH + CH11 + CHe 

rE ¥ ch G7. =2,) GD TO 70 
HTTH = @ S#¥CSTHACTH)D 

CAT = CLeHTTH*PS1#6011 

CIT = CEeHTTHePSe4#Oe 

CY = CY © CLT © UST L-w 


a ee IR 


rT 
" 


5 u L 
CWE 2 Ete CVI t+CGeFTLI*PS1 - €C4 ) 
CHE = Cts HE=-COPPTISePSS = ECA 2 
CALL FTE (S2 FT 2 
c C34 WES+CG@RPTSO FSS - COZ 
; $4 t-CO*FT2)*PS2 - CCF > 
FT eCGHEFT SHES 2 
ALL PT TS4, FT4 
eT = 4 CMAP Ze Peg 
CY = CY +#¢CV4 + CVE - CEFT ~ C4FT) #04 
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FUNCTION C(Z) 


CZ) = O.5*(1 + ERFIZ)) 


* ERF(Z} IS EVALUATED USING A RATIONAL PCLYNOMIAL APPROXIMATION 


* REFERENCE (HANDBK MATH FUNCT BY ARRAMGWITZ ANC 


7 SECTION 7.1226) 


IF ( Z.GT. 22) GE TO 10 
IF ( ZeLTe-2-) GC TO 20 


AZ = ABS(2) 


a) 
" 


1L.0/(1.0 * o4704T*AZ) 


Y = 160 = P#(63480242 — PR(.O958TYE 


IF (2) 29496 


© = (1.0 ~ Y)/2. 


nN 


RETURN 
40 = .5 
RETURN 
60 = (1.0 + Y)/2. 
RETURN 
16 © = le 
RETURN 
20 G = Of 
RETURN 
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L.7 STEPPED - CYLINDER 


The far-field scattering from a stepped-cylinder target 
configuration shown in Figure L.4+l has been formulated using 
the basic Ruck-Ufimtsev technique and a rough approximation 
of the multiple reflections of the field components between the 
cylindrical and flat surfaces (Ref. 8).. The basic scattering 
formulation is the following: 


e(8)ivi= Fym { g(11) + (12) + g(1) + e(2) + 8(3) + B(4) 
+ 8(9) + 8(5) + 8(6) + 8(7) + 8(8) + (10) 
where g(m) represents the sum of the fringe wave scattering 
and physical optics response associated with 
edge m. 


For 0<@<7/2, 


g(12) = PUL { 33444 [c(1.5) ¥ B(1.5,7/2-0)] - C(1.5)JJ93} 
g(1l) = im LID; [c(1.5) ¥ B(1.5,7/2+8)] - c(1.5)JJo,} 
g(2) = aye P? (ay) [c(a.5) F-B(1.5, 9/248)] - c(1.5)I5,, | 
g(4) = tdi [cc1.5) ¥ B(1.5,7/2+0)] - c(2. 5)JJy,} 
g(9) = age Lie. [cc.5) ¥ B(2.5,0)] - c(1.5)JIIq4} 
g(l) =+a = 0.5 JJ,,. [tano + cote] } 

g(3) =+a,e [0,5 33,5. [ tane + cord] | 


; a 
g(6) = a,e P2 Rit, [cca.5) F B(1.5, 7/2-8)| 


. C(1.5)I3 99 | R61, 
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Fig. L.7-1 GEOMETRY OF STEPPED=CYLINDER 
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g(8) = ye Ps {23 [c(1.5) ¥ B(1.5,7/2-0) 


‘ eG 7a, QgIQs 


ip, 
g(5) = 4;.¢ { 0.5 33454008 | Q, 19, 
ae ip 
g(7) =+ age 4 | 0.5 JT, g C080 } aside 
For 7/2<@<7m7 
g(10) = a,e 2 Matures [cc1.5) ¥ B(1.5,0-7/2)| - c(1.5)34,, | 
g(9) = ae ee [cq. 5) + B(1.5,22-6) - C(1.5)IIz4 } 


+1 


ae 
= 4 
g(4) a,e (33, [c(1.5) 


F B(1.5,7-0)] - c(1.5)JJ,, } a 


“4 
ip a 
g(2) = a,e 2 t 3 5 [c(1.5) + B(1.5,7-0)| P C(1.5)JJ5, | Q, 


; hae 
g(1l) = aje Fil { 334). [cc1.5) ¥ B(1.5,7-9)] - c(1.5)I55; QQ 


i 
g(1) =+a,e *f 0.5 JJ tano} Q 


1l- x 


g(3) 


where the upper and lower signs in the previous equations 
correspond to the case of vertical and horizontal polariza- 
tion, respectively, and 


sin = 
C(n) = aa I—— | 


1 
cos ~-l 
n 


" 
1+ 


ip 
aye 4 { 0.5 JJj9. tand } Q3 


sin n 1 
sds. ee foal ae 


cos — =- cos —= 
n n 


JJ =: ™ J,(2ka_sin@) + J, (2ka,sin@) 


JJon = J, (2ka,sin@) + Jo (2ka sin@) 


oe eae 2k (hy + ho) cos®@ 
Pp = - 2k hy cos@ 

a 2k hy cos® 

Ps = 2k(hy + hy) cos 


Q ig) 7 Q¢2ka 4, Cay, = @)) 5 Qy = Q(2ka, (F+8-0)) 


(4) (3) | 


(22) g= eles sec ds ; Q3 = Q(2kaz(F+8-8)) 


2 2 2 7 
§ = D 
Q, = Q(2ka, (-8)) 
Qn = 1 for Q,, 70-5 Sy 
0 otherwise 
a,-a 
-l/ 2 1 
@, = tan ( h ) 
l 
ao am. 
a, = tan’ ( 2hp ) 
a,-a 
- 4 1 
a, = tan"}( 
3 h7zh,) 


Q(z) = 0.5 [1 + erf(z) | 
a =a, + hy tand 


a =|a, + 2h, tan@, for @< ay 


2 


a, + (2h, +h, ) tan@,for a,< <a, 
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Most of the equations can be obtained in a logical and 
straight-forward manner from the well defined procedures of 
the Ruck-Ufimtsev high-frequency scattering technique. Un- 
fortunately, no first-order technique can be used to obtain 
an accurate formulation of the scattering from surfaces and 
edges which are partially shadowed or just within the shadow 
region. The screening functions of the Ruck-Ufimtsev tech- 
nique have been used in computing the response of several 
aerospace targets at aspect angles in which some of the 
target surfaces and edges are shadowed; however, the use of 
these functions is a first-order approximation, and multiple 
diffraction and reflection must be considered in order to 
accurately describe scattered returns from shadowed target 
regions, 


The formulation of the multiple scattering between the 
surfaces adjacent to the concave edge structures is expressed 
in the following equations , where the "*'' is used to distinguish 
these expressions from the, first order approximations on page L-126: 
* 4 
g(1*) = Vak ay e Pl 


~ 


* 
V4k A oP3 


g(3*) = 
* %* ipl; 
g(13) = V4k Al3 e 
where Ay h, sin@ 
= ja sin 
l ry for Q<a, 
= V a] (ag-a,) cos® e2 Qa, 
* > 
Ag = Va, 2h,sin@ for Q@<a, 
= Va, (a,-a,)cos® Q2a, 
k= h i - = < 
AL; va,{ 1Sind-(a,-a, )coso] a, S0sa, 


= Va, {(a,,-a,)cos-2h, sing} Ge <@< a, 
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* 


es 3 a ee ern nn | 


‘* 
and ee -2k(a, sind + h,cos@) 


3 = -2k(a,sin@ - h,cos®) 


2 


* 
P13 * ~2k(a,sin® - hycos®) 

This approximate representation of the multiple reflected 
fields cannot be considered to be an exact or rigorous formu- 
lation of the scattering mechanisms involved. Nontheless, 
this preliminary analytical expression is simple, and the 
computed results do provide a realistic measure of the 
target scattering. 

L.%. i ‘Inputs 


The subroutine inputs are read from cards or passed in 
common blocks. The parameters passed in unlabeled common 
include: 

NMIN = minimum frequency sample number 

NMAX = maximum frequency sample number 

DF = frequency increment (in MHz) 


FC = carrier frequency or center frequency (in GHz) 


The parameters passed in a labeled common block 
include: 


ASPECT = azimuth angle (in degrees) 


ITT Read data option set to 1 or 2 
= 1 Read target dimensions from input card 
= 2 Use dimensions input on prior read 


The card inputs are the following: 


a* 


MATH 
SYMBOL VARIABLE DEFINITION COLUMNS 
Card L ay Al Smallest cylinder 1 - 8 
radius (inches) 
ay A2 Middle cylinder 9 - 16 
radius (inches) 
ay, A4 Largest cylinder 17 - 24 
radius (inches) 
hy HL Length of smallest 25 - 32 
cylinder (inches) 
hy H2 Half-length ‘of middle 33 - 40 
ie cylinder (inches) 
hg H3 Length of largest 41 - 48 


cylinder (inches) 
L.7.2 Outputs 


The data base output consists of four linear arrays, 
EVVR, EVVIT, EHHR, EHHI, which contain the real and imaginary 
parts of the vertically and horizontally polarized back- 
scattered fields (in meters) at frequency increments of DF 
MHz from NMIN*DF to NMAX*DF. 


L.7.3 Restrictions 


L.7.3.1 Physical Dimensions 


All target dimensions should be large with respect to 
the wavelength of the illuminating field. In formulating 
the Q functions to approximate the effects of ergo on 
an edge or surface, the assumption that @ 1<% was 
utilized, since this angle relationship was Bibi ed in the 
geometry of the target for which the formulation was devel- 
oped. If this angular relationship is not present, the 
arguments of the Q functions must be modified. 


L.7.3.2 Output 


The output arrays are passed in the argument list and a 
value is computed and stored only in locations NMIN to NMAX. 
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L.7.3.3 Restrictions 


The azimuth angle is restricted to the region between 
0 and 180 degrees. In addition, the specular azimuth of 
0, 90, and 180 should not be used. 


In order to compute the response of these angles, an 
angular offset of approximately 0.01 degrees should be 
used. 


L.7.4 Definitions of Selected 
Terms Used in Subroutine 


SV12 = C(1.5) # B(1.5, 7/2-0) for @ <7/2 
1p = 
PHIP = JJ,,, 
= JS) = J, (2ka_ sind) + J, (2ka_sin®) 


where m = 1 and the lower (+) sign is used 
BCIX = r 
C1IxX c(l 3)JJ55 


where BCl = JJ54 


= JJ, = Jo(2ka,sin®) + Jo (2ka,sin@) 
where m = 1 
cll = e4P11 
where 


Pit = - 2k(hy + hy) cos®@ 
i _ 
CV12 = aje P11 {53114 [c(1.5) ¥ B(L.5, 72-0)] = C(1.5)JJQ}} 


for @ < 7/2 


ove = g(6) = aye 2 { sy, [C(L.5) ¥ BCLS, 7/2-0)] 


-C(1.5)JJy, } 619, 


| PSDL = p, * -2k(a,sin@) + h,cos@) 
ADL = Ay * Vajh sin® for @< &,,and 
ADIL = Ay = Va, (a,-a,)cos® for Q> @, 
L.7.5 Subroutines Used 
Subfunctions: 


1. Q(X) returns the value of the exponential smoothing 
function. 


Subroutines: 


e l. BESL(XCx, BJxQ BJxl, BJx2) computes the Bessel functions 
of order 1, 2, and 3 for real argument XCX and returns 


Jo (XCx) in BJx0 
Jy (XCx) in BJxl 


J5 (XCx) in BJx2 
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EHHI 3 


(MODEL 435 


SEY 


FC. 


TONS 


NS of 


H1, 


TECHMI GUE 


Fill. 


ia] 


NOT 


He, 


Ta 


Hs 


PARAMETERS’. “07. % 


“yE?. €a7% 
4, 


” 
‘ 


4 


nl a 


' 
fi 


) 


aA 
+: 4: 


CHANGE 


ChHd, CV 
He. Geo 


AD-A031 440 GENERAL DYNAMICS FORT WORTH TEX CONVAIR AEROSPACE DIV F/G 17/9 
ENDO ATMOSPHERIC EXO ATMOSPHERIC RADAR MODELING. RADAR CROSS SE=-ETC(U) 
JUN 76 R J HANCOCKs F H CLEVELAND eueene ae 

UNCLASSIFIED RADC-TR=76-186-VOL-4-PT=2 


/ ALFS ATANCF44“THEPH1) 


C 
FiNAtLe2 = PIi<Alre 
PIMALS = PI=ALF2 
c 
810 = ALFL*kTO 
Seb = ALFEerRToO 
NEG = ALPERT 
C 


WRITE (6, €GQ5? 840, 82D, M30 
6G05 FORMAT | ALFA = ’,F8.3,/¢2" ALF2 © ’, FS. 3.7.’ ALES = °, FS, 2 


TKFC = FITK#*FC#e 
SN S eS 
KNOe = ®N'e 
SOxN = SINCPICSNICAN 
EXN = = SON N 
CRON = COS¢PICMN) 
E 
& CONTINUE 
TH = ASPECT#DTR 
TH = SINCTH) 
TH = £LOSCTHS 
TANTH = STH/£LTH 
CK! = @ S#TANTH 
CKY = @ S/TANTH 
CLC! =x OR LECEY 
iF ¢ TR Gi PIQ2> oh To 26 
Bic = SORN-CCFEON ~- COSCCPIOG-THI “ENGe) dD 
Gi4 = SORNA CRON - COSCCRIGe+THI RNG) D 
BS = SOMNCOCPON - COSCTHARNOZS 5 
EOGE DIFFRACTION COEFFICIENTS ¢THETA. LT. PIl/¢2) 
Suge = CKN- Ble 
SHie = CAN+EL 
S¥Vil = CHN- E11 
SH14 = CHN+B11 
ovo = CRN-EBS 
SHS = CHN+ES 


C EFFECTIVE AREA COCUBLY REFLECTED > TERMS CALS? 
TF CTH=ALELD Si. 52,52 
SL ACL = SAL*HL*STH 
GO TO Sz 
Se ACL = SAL#A21*CTH 
= CONTINUE 


IF ¢TH-ALFe? 84,55, 8 
"4 Abe = SAZtTHZ*STH 
Gc Ta 56 
SS Abe = SAe*AdZe*CTH 
S6& CONTINUE 


an 


IF <TH-ALF4A) S77, S58, 5 
S& IF <TH-ALFe) SY, 59, 57 
S93 IF «TH-ALF2) 64, 64, 64 


L-/46 


ol 


ADS = SAL®CHL*STH - A&L*CTH) 
Go TO 6S 

ADS = SAL*#CAde# CTH - THeaSTH> 
Go Ta 63 

ADS = @ 

CONTINUE 


ASSUMED THAT ALFS .GT. ¢ALF1 OR ALF] D 


SMOGTHING FUNCTIONS AND ILLUMINATE AREA USING SLIDING SHADOW 
BOUNDARIES : 


ISWAZ = 4 
IF < TH GT. ALFE > GO TO 2s 
ISWAS = 2 

06 9 = MK TKEORAZHCAILFA-TH)D 


IF (06 LE. @&5) G6 = @ 
te = OCTRFO#A4+€ALFE-TH) > 


Ir €@o . bE 6.8 9 GS. = -& 

ASS = Ad + H1*TANTH 

IF «< TH . GT. ALFA > GO TO 14 

AS? = Ae + THe TANTH 

Go Ta 2& 

AS? = AL + THEPHL+TANTH 

GOT 25 
THETA .GT. FI/2 
B1G@ = SORN/CCRON - COSCCTH-PIGZ) “¥NGZ>) 
ES = SOANSCCRON - COSC CTRPIG2-TH) “Ande ) > 
G14, = S0AN/CCRON - COSCCPI-THI CHENG) D 
<DPGE OLFFRACTION COEFFICIENTS ¢THETA. GT. FI/¢2) 
S¥1ia@ = CKN-614 
SHia@ = CHN+E414 
S¥o = CXN-EBS 
SHS = CHN+ES 
SV4L CeN-EBAL 


SH141 = CKN+E44 

WALE = PIOZ+OLT-TH 

SMOOTHING FUNCTIONS 

GL = G@CTKEC#AL*eX ALP) 

2 = Ge TKE CASH CP IMALZ-TH))D 
C2 = GCTKFC#AZAKALF) 

C4 = OCTKFCAAd#CPI-THDD 

G11 = OC TKFCHAL*CPIMALE-THDD 
CONTINUE 


FREQUENCY LOOF 


CO 36 I = NMIN, NMA 


KI = I-14 

AK = XIeFTKDF 
AK2 = AK+AK 
AK2S = AK2*STH 
KCL AL*AKZS 


= 
KCZ = AZKAKZS 
= A4KAKZS ust) 


mmr 


mw | 


GENE F 


FRASE 
POL 


COUEL 


CAT 


oe 


BESL¢XCL, BU1G, BJ11, BJ12) 
BESL¢XCe, BJ2@, Ble1, B22) 
BESL¢XC4, BJ40, B41, BI4e) 


= BILG+bI12 

= BIZa+EJe2 

= BI4G+b Ide 

= BCL#*CKN 

= EBCL#CRN 

= EC4#CKN 

= CMPLACEILG, BY113 

= CMFPLACEI26, BY 21) 

= CMPLECET4@, B41) 

= CONJGCFHLF) 

= CONIJGCFHEF) 

= CONIGCFH4F) 

= ~-AK Se CHe+H1o4#CTH 

= -AKS#He*OTH 

= Ak ee CHS+HEo*CTH 

= CMPLECCOSChS129, SINC PS44 95 
= CMPLACCOSCRSS3, SINCPSS) 
= CUNIG¢@ Ce 

= CMPLMCEOUSEPSSS, SING PES 73 
TH -@fo “tae + Gh Tess 

AL REGION <THETA .GT. FIL FPS oFNG. THETA 2ELT. FI¢2) 


OF DUUBL + REFLECTED SCATTERING TEENS 
= -AK2*CALesTH + Her lTH) 
= -RAike2eCAe#STH - Hea CTH? 
HoreTAy 


»SING ESOL) 3 
2s SINCPSOS) 3 
i, SINC ESOE) 3S 


Y REFLECTED SCATTERING ¢TOTAL? 
= SGRTCAKS+AK 22 *CADL+CO1 + ADS*Cbe + ADS*COS) 


UPIMTSEY SCATTERING (THETA. LT. $@,E0GES ILLUMINATED AT 9@) 
ALCL i eo S¥ie4+PHIP-BC1K) 

ALCL ¢SHL24#PHIF-BO1k) 
AL*C11#¢SVY41*PH1-EC1#) 
AL*C11i#¢SHL1*FHI-BCI x) 

ASe Cee CSVi1ePHe-EC2x) 

AeeCeecSHi1ePHe-BC2]N) 

Adeldee SV14#PH4-BO4 ND 

Ad#C de (SHIT #PH4-BC4) 

Ad*eCSec SYGRPH4-BC4e ) 

Ad#eC Sa SHOePH4-EC4 

Al+*eCowPHi+C1 Poy 

AZeC4ePHerC1Pcy 

C¥YLEF+CV114CV2+CV Geto +04 +0VE L-/39 
CHL2+CH1144+CH2+CH44+CHS -CY41 -CVE 


Be QR ae eau es 


— 


40 


wt 
2 


Sy = SY + CADSFI 

SH = SH + CADS/SFI 

GO TO (46,34), ISWAZ 

CGR TINUE 

ACINUTH NERF ZERG, ACD IN FARR EOGE TERMS WHICH 


DIRECTLY ILLUMINATED 
RCSS = ASSeAK es 
RCS? = ASP*AK SS 


CALL BESL(XCSS, 656.651, &82) 

CALL BESL¢(XCS?, 676,671, B72) 

PHSP = CMPLXCES@, B51) 

PHFP = CMPLXCE7@, B71) 

CVE = AeelCeecPHZPeSy1 o-BCoN) #0 
CHE ALS Cle PHEPeSH1e-EBCox > #6 
CVE = Ad#C4eCPH4PeSv1e-BC 4x) #08 
CHS = Ad4#+C4eCPH4Ph*eSHi12-EC4 ko #08 
CYS =-ASS#CEePHSPe Ch vee 

CYP =-ASPHeRO4HPHF PECK VY #OS 

sy = CVE + CVE +CVE4+CVF7 + SY 
SH = CHE + CHE -CYS-CY? + SH 
Go Ta 46 

THETA .GT. FIre2 

CONTINUE 


RUCK-UFIMTSEY SCATTERING (THETA. GT. 9@3 


C¥Y1@ = A4*CoeCPHASP#SV1A-BCdx) 

CH1G = AdtOC Ge C PH4+P#OH1G@-EC4 > 

CVS = Ade CoeCPH4eSyO-EC4 kD 

CHS = Ad#Co#eCPH4 tOHo-BCd x) 

CV¥4 = Ade Ge CPHdeSV14-BO4N) #04 

CH4 A4#C+eCPH4tSH11-BC 49 eld 

CVS = AS*CSkcPHS*SV11-BC2k) #02 

CHE = AdeCeeCPH2#SH11-BCex #2 

CV1i1 = AL*C11* CPHL SV11-BC18 #014 
CHL = ALeCLi eo PHL eSH1L1-BCLR #041 

CVi = ALeCe#FHi CK i 04 

CVT = AgeC4ePHee CK 1 OF 

cy = CYLG+CVS4+CV G+ CVeeCV4d +0Vi+CVE 
SH = CHLG+CHS+CH4+CHe+CH11 -CV1I-CVE 


CUNTINUE 


-SYaSPIE 
SH*SPIK 


' 
sy 
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cH = 


REFERENCE FHASE TC FRONT OF SECOND CYLINDER 


V 


S 
SH 


SVECS 
SH*EC4 


EVVRCI) = REALCSY) 

EWVICI> =-AIMAGCSY> 
EHHR¢I> = REAL ¢SH> 

EHHI¢I) =-AIMAG¢( SH) 
CONTINUE 

RETURN 


L-\34 


ARE 


as aati a rig gneyeeeeeeee 


ENG 

FUNCTION @¢2) 

2> = & Seth + ERF CZ) 

+ ERFOD) IS EVALUATED USING A RATIONAL POLYNOMIAL APPROXIMATION 
* PEFERENCE ¢HANDEK NATH FUNCT BY ABRAMOWITZ AND STEGUN, 

+ SECTION 7. 4. 26) 


2.) GO TO 16 
2.2 GO TG 2e 


Dn m~ 
mA 
mrynm 


1@ + . 47@d7*Ae>d 


C B4ehe¢e - Pee B98G79S - F447 ESSEeP DD eKEXP CO -AZ* AZ) 


on 


422 .% 
PETURM 
60 = 1.84 © VI 
RETURN 
1@ @ = 1. 
RETURN 
“6 0 = @ 
FETUFN 
ENT: 
SUBRGUTINE BESL ¢ &, E&, BI, Be 3d 


+ BESSEL FUNCTICN SUBROUTINE UTILIZING POLYNOMIAL APPROXIMATIONS 

+ COMPUTES J@. 04, 0R J2 FOR POSITIVE REAL ARGUMENTS 

+ REFERENCE CHNDEK MATH FUNCT Bt ABRAIMGWITZ AND STEGUN SECTION 9.4 ) 
S= 41.6 
IF <x LT. @ @> S=-4.6 


= ABSCK? 


EF C8 2 GF. Sbeey Gd Ee 3 


B@ = 41.@ 
Bi = @. 6 
bo = & & 
| sR SS 
RETURN 


© CONTINUE 


i 


¢ * .8& 3.9 86 Te s 


B= 1 + NAHC-2. 24999974 HLL. ZESEZOE+ KC -. SLEFEEE+ NA*C. 444479 


i+ MLe¢-. @@59444+ Ki*e 16-4 9))0) 2D 
60 TO 16 
9 KZ = 3.7% 
Fa = TUTES4SE +h 2HC-. PTE-E +KeeC-. GOSSS7AG +KSHC-. GSLZE-4 +K2% 
1 (. @B137237? +K2eC- PAEOSE-F +K24G. 14476E-3F 999) OD 
TH = K - | 7PESFIGLE +K2HC-. G41GETG7T +X2eC-. TOS4E-4 +KAeC. GHSEZS73 
1 there ~ GOBS4125 +X2eC- GOG29TTT +K2*O BOGLZESE 9d) 


B = Fe*COS¢ TB) /SORTCX) 
hi Ra = & L\40 


c 
21F (x .GE. 3. > GO TO 15 
Xi = X73 
Ki = Kiexd 
B= KeC 5 +X1e(-, 56249985 +X1*(. 21093573 +K1H(-. @2954289 +K16 
1 (.@G443319 +X1#(-. Z1764E-3 +X1#@. 1109E-4)9))) > 
Go TO 26 
c 
19 K2 = 3. /x 
Fi = . 79788456 +X2«(. 156E-S +XZHC. O1659667 +KZHC. GOOL17105 +K2e 
1 (=, @8249541 +K2*(. OG4113653 -. OBOZOATT*X2 9)))d 
TL = X - 2. 35619449 +X2e(. 12499612 +XZ2e(. SESE-4 +K2e(-. 0627879 
1 —»- #K2#C. @BE74I4S +X2H(. BOE7IEZ4 -@. BOAZS1E6*KZ >) 
B = F1ieCOS(T1)/SORTCX) 
c 
20 B1 = B® S 
K =X # S 
B2= (2. /X)*B41 - BO 
56 RETURN 
END L~\4\ 
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